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Abstract 

We systematically calculate the relativistic corrections to the polarization variables of prompt 
J/if photoproduction and hadroproduction using the factorization formalism of nonrelativistic 
QCD. Specifically, we include the and 3 P^ color-singlet and the 3 S'j 8 ^, an d 3 P® color- 

octet channels as well as the effects due to the mixing between the and 3 .D® channels. We 
provide all the squared hard-scattering amplitudes in analytic form. Assuming the nonrelativistic- 
QCD long-distance matrix elements to satisfy the velocity scaling rules, we find the relativistic 
corrections to be appreciable, especially at small transverse momentum px and large inelasticity 
z. The results obtained here and in our previous work on the unpolarized yield [Phys. Rev. D 90 , 
014045 (2014)] will help to render global analyses of prompt J/if production data more complete 
and hopefully to shed light on the J/if polarization puzzle. 

PACS numbers: 12.38.Bx, 12.39.St, 13.85.Ni, 14.40.Pq 
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I. INTRODUCTION 


Through concerted efforts from both the experimental and theoretical sides, we have 


deepened our understanding o 
years (for a review, see Ref. 


' the mechanism of J/-0 production at high energies in recent 
1] and references cited therein). Yet, there are still major 
challenges to the theoretical models. One of them is the long-standing J/^j polarization 
puzzle. The polarization of promptly produced J/ij) mesons, which are produced either 
directly or via the feed down from higher excited states such as XcJ (J = 0, 1 , 2) and i/i' 
mesons, was measured by several experiments in different environments, including CDF at 


the Fermilab Tevtron 


at BNL RHIC fl; and ALICE J8j, CMS 


, 3 ; 


HERA-B 


4], ZEUS Q, and HI 3 


and_Hl [6| at DESY HERA; PHENIX 
at CERN LHC. Unfortunately, 


9], and LHCb |nj 

these measurements could not yet be explained by theoretical analyses in a way consistent 
with the world data on the unpolarized J/i[) yield. The production of heavy-quarkonium 
states, like the J/-0 meson, involves the creation of the heavy-quark pair {QQ) and its 
subsequent transition to the hadronic bound state. Different ways of treating these two 
parts result in different models. Polarization variables are more sensitive to the fine details 
of the hadronization of the QQ pair than production rates and, therefore, play a key role in 
testing the theoretical models. 

Due to the fact that the heavy-quark mass rriQ is much larger than the asymptotic scale 
parameter of QCD Aqcd, the heavy-quarkonium state can be approximately treated as 
a nonrelativistic system. The factorization formalism based on the effective field theory 
of nonrelativistic QCD (NRQCD) llHT3j is nowadays considered to be the most favor¬ 
able theoretical approach to study heavy-quarkonium production and decay. In this frame¬ 
work, the heavy-quarkonium production cross sections are factorized into process-dependent 
short-distance coefficients (SDCs) and supposedly universal long-distance matrix elements 
(LDMEs) |l3. The SDCs, which describe the production of the QQ pair at energy scales 
2 rriQ or larger, can be calculated perturbatively through expansions in the strong-coupling 
constant a s . The LDMEs, which measure the probability of the hadronization of the QQ 
pair, are weighted by definite powers of the relative velocity v of the heavy quarks in the 
heavy-meson rest frame l^. In this way, theoretical calculations are organized as double 
expansions in a s and v. Quantum corrections come as terms of higher orders in a s and rela¬ 
tivistic corrections as terms of higher orders in v 2 « ot s . Thus, next-to-leading-order (NLO) 
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results include corrections of orders 0(a s ) and C?(v 2 ) relative to the leading-order (LO) 
result. Next-to-next-to-leading-order results also include terms of relative orders 0(a 2 ), 
0(a s v 2 ) , 0(v 4 ), and so on. A crucial difference between the NRQCD factorization formal¬ 
ism and the conventional color-singlet (CS) model is that the former allows for the QQ pair 
to be in any possible Fock state n = 2S+1 L ^ with total spin S, orbital angular momentum 
L , total angular momentum J, and color multiplicity a = 1, 8, while the latter is restricted 
to the CS state, with a — 1, sharing the quantum numbers S, L , and J with the heavy 
meson. The color-octet (CO) mechanism, i.e., the appearance of CO states, with a = 8, is 
a distinctive feature of NRQCD factorization. 

Two decades after the introduction of NRQCD factorization jl^], all the relevant observ¬ 
ables of prompt J/ij) production are available at NLO in a £ . In particular, these inc 


yield 


14] and polarization 15] in photoproduction, yield 16-tl8| and polarization 


in hadropro duct ion, and observables in other production modes 


24 


inc 

ude 

19 

-23] 


26]. Since the LDMEs 


of the T) c meson are related to those of the J/-0 meson via heavy-quark spin symmetry, the 
Tj c yield 27-[29] must be included in this list, too. The J/-0 LDMEs determined at NLO in 
NRQCD via a global fit to measurements of the unpolarized J/V’ yields in photoproduction, 


3a 


20 


21 


3 . 


hadroproduction, and other production modes 
in hadropro c 
the LHC |9 
data 


26] lead to predictions of J/-0 polarization 


uction |l9j that are incompatible with measurements at the Tevatron 2, 3] and 
On the other hand, J/ip LDMEs determined only from hadropro duct ion 


are incompatible with data from photoproduction and other production modes 


3d ]. In other words, the 


ong-standing J/ijj polarization crisis of NRQCD, which had already 
been observed at LO 311, is substantiated at NLO in a s . This crisis has recently been ag- 

n nnnn 

gravated by the observation [27j that all the up-to-date J/^ LDME sets |20|, 21|, |26|, |32] lead 


to NLO NRQCD predictions for the r/ c yield in hadroproduction 


the recent LHCb measurement [33]. LInfortunately, attempts 28] to reconcile the measured 


hat significantly overshoot 


y c and J/ij) yields of hadroproduction by resorting, to pre-LHC LDMEs 17] fail to describe 
the J/ij} polarization measured at the Tevatron 23] and the LHC 28] in the regime of large 
Pt values and central rapidities y. The analysis of Ref. j^ij also misses its goal of achieving a 
coherent description of the J/-0 and yields and the J/ijj polarization in prompt hadropro¬ 
duction. Obviously, NRQCD factorization presently faces severe challenges with regard to 
the predicted universality of the LDMEs at NLO in a s . 

To resolve the J/-0 polarization puzzle and to clarify the universality problem of the 
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NRQCD LDMEs, we systematically calculated, in our previous work 34], the relativistic 
corrections, of relative order 0 (v 2 ), to the unpolarized J/0 yields in photoproduction and 
hadropro duct ion by including both the direct and feed-down contributions. We found that 
the 0(v 2 ) corrections are appreciable, except for the CS channel of hadroproduction, and 
that their line shapes greatly differ between photoproduction and hadroproduction, which 
may offer a chance to improve the goodness of the state-of-the-art determinations of the 
LDMEs. In this paper, we take the next logical step by calculating the 0(v 2 ) corrections to 


the J/-0 polarization observables in photoproduction and hadroproduction. A' 


14 


1 these 0 (v 2 


23] to render 


corrections must be included on top of the respective 0 (a s ) corrections 
the NLO predictions complete. 

We organize the remainder of this paper as follows. In Sec. HH we shall briefly describe 
how to calculate the SDCs for the polarization parameters in the direct and feed-down 
processes of photoproduction and hadroproduction. In Sec. m we shall present and discuss 
our numerical results. A brief summary will be given in Sec. IIV1 All the relevant analytic 
expressions will be collected in the Appendix. 


II. NRQCD FACTORIZATION FORMULA 


The polarization of the J/0 meson is measured through the angular distribution of its 
leptonic decay, J/ij) —> l + + l~, which may be parametrized as 


W(9, 0) = 1 + \ e cos 2 9 + A<£ sin 2 9 cos(20) + A 0 $ sin(2 9) cos 0, 


( 1 ) 


where 9 and 0 are the polar and azimuthal angles of lepton l + in the J/0 rest frame, 
respectively, which depend on the choice of coordinate frame. For example, in the recoil or 
s-channel helicity frame, the z axis is chosen to be the J/0 flight direction. The parameters 
\q, \<p, and A 04 , in Eq. (pQ) may be related to the helicity density matrix dc r A y, which describes 
the interferences between the states of helicity A and X' in J/0 production 35] 


A fl — 


do 1,1 — dcro.o 


A ^ — 


do i,i 


A a a. — 


\/2 Re d(j\ 


( 2 ) 


d(7 11 T dcr 0i0 d(Ti \ T do 0 0 doi ^ T do^ o 

Invoking the Weizsacker-Williams approximation and the factorization theorems of the 
QCD parton model and NRQCD 131], we may write the hadronic helicity density matrices 
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of both the direct and feed-down processes in the general form 


daw (A + B —y 



dx 1 dy 1 dx 2 fi/A{xi)fk/i{yi)f j / B {x 2 ) 


x da xy (k + j H + X)Bt(H ->• J/V> + X), 


( 3 ) 


where f t /A (x) is the parton density function (PDF) of the parton i in the hadron A = p,p or 
the flux function of the photon i — 7 in the charged lepton A = e~, e + , fj/i(y 1 ) is Sij5(l — yi) 
or the PDF of the parton j in the resolved photon i — 7 , Bt(H —> J/'ijj + X) is the branching 
ratio for H = J/'0, XcJ, '0 / with the understanding that Br (J/ip J/^ + X) = 1 , and 
ddw (i + j H + X) is the partonic hclicity density matrix. Through relative order 0(v 2 ), 
the latter may be factorized as: a 


dd X y(i + j H + X) = 


n 


dFt\,(n) ,^tt, .. dG l A,(n) .. 

iA y o" n + A - (P (n)) 




l O(n)~ 


m { 


l V(n)~ 


( 4 ) 


where O h ( n) is the four-quark operator pertaining to the transition n —> H at LO, with 
dimension do(n)i V H (n) is related to its 0 (v 2 ) correction and carries dimension d-p pq = 
do(n) + 2, and F^ x ,(n) and G xx ,(n) are the appropriate SDCs of the partonic subprocess 
i + j —> cc(n ) + X. The definitions of the O and V operators and some of their properties 
may be found in Refs. jl3, 34] . Working in the fixed-flavor-number scheme, the parton i 
runs over the gluon g and the light quarks q = u,d,s and antiquarks q. 

As in the unpolarized case [34j, the relevant partonic subprocesses include 


9+7 -> cc( 3 s! 1 A V 4 Vb 8 | )+ 9 , 
9(9)+7 -> cc( 3 Sf 1 , 1 Sj 8| , 3 P| il ) + 9 ( 9 ). 

9 + 9 -> cc( 3 Sl 1 ’ 81 , 1 Sf 1 , 3 P] 1 ’ 81 )+g, 
9®+9 -> cc( 3 Sf l , 1 Sj, S| , 3 Pj 1,8| ) + 9(9). 
9 + 9 -+ cc( 3 s| 8| , 1 4 8| , 3 pl 1 ' 8| )+9. 


( 5 ) 


The SDCs F xx , (n) and G l xx ,(n ) may still be calculated with the help of the spinor projec¬ 
tion method as explained in Ref. 34] when the polarization four-vector e M (A) of the cc pair 
is kept generic. In the following, we only explain the differences in computing the hclicity 


The spin-flip interactions are 0(v 3 ) suppressed 36|, so that NRQCD factorization still holds for the 


helicity density matrix through relative order 0(v 2 ). 
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density matrix elements with respect to Ref. 34], We only need to discuss the nontriv¬ 
ial cases of S — 1. In the case of direct J/ij) production, the total spin S of the cc pair 
can safely be identified with the spin S of the J/-0 meson through 0(v 2 ) because spin-flip 
effects occur only at 0(v 3 ) 36], and the index L z related to the orbital angular momen¬ 
tum of the cc pair may be summed over. Starting from the partonic scattering amplitude 
M\(n) = M(i + j —> cc(n, A) + X ), where A is the helicity related to the total spin S of the 
cc pair, the helcity density matrix element is thus evaluated as p xx ,(n) = J2M\(n)M x ,(n) by 
summing over the orbital angular momentum L z of the cc pair and the spins and colors of 
the other outgoing partons, contained in system X , and averaging over the spins and colors 
of the incoming partons i and j. Through 0(v 2 ), the general Lorentz structure of p l xx ,(n) 
implies the decomposition 


P\\'( n ) = l A k ( n ) + ^ B k ( n )] «T<(A)e„(A0, 

k 


( 6 ) 


where k labels the process-independent four-tensors 

J lIJ — — U l b 1 ' J llJ — J w — j _ k ,l k v 

*1 ~ > *2 '■'l '■'l j *3 — ^ 2 ^ 21 *4 — M ^2 ' ft '2 ft, lJ 


( 7 ) 


with k\ and A; 2 being the four-momenta of the incoming partons; A l £ (n) and B^ (n) are the 
SDCs, which are functions of the partonic Mandelstam variables s, t, and u and depend on 
the considered partonic subprocess, but not on the choice of coordinate frame; and q is the 
relative three-momentum within the cc pair. Because P ■ e(A) = 0, where P is the total 
four-momentum of the cc pair, and P 2 = 4 E 2 34], e M (A) also implicitly depends on q 2 . To 
obtain the complete results at 0(v 2 ), the contractions between and e*(A)e„(A') also need 
to be expanded as series in q 2 after choosing a coordinate frame. Writing 


^;(A)b(A') = 4, + q 2 4 + 0(q 4 ), 


( 8 ) 


it is straightforward to obtain 


EMU 

dotn )—4 

m c 


1 

2s 


dLIPSj]4«(n)< 


AA') 


= A f dLIPS V [4>) (Kc\ y + d\ y ) + B«(n)c‘ v ] , 


m c — - k 

where the definition of the factor K is the same as in Ref. 
frame dependent. 


( 9 ) 


34]. Note that c xx , and d xx , are 
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The radiative decays \ c j —>• J + 7 are not only affected by the El transition, but also 
by the higher-multipole M2 (xd and Xci ) and E3 (xa) ones. A detailed study Q revealed 
that only the angular distribution of the photon is very sensitive to the higher-multipole 
contributions, while the latter are negligible for the angular distribution of the subsequent 
J/ijj —y l + + l~ decay. Therefore, it is sufficient to calculate the helicity density matrix 
elements of J/'ijj production from \ c j feed down by connecting the J/^ polarization four- 
vector to the total spin of the cc pair in the helicity density matrix elements of XcJ production 
and multiplying the outcome by the branching ratio of \ c j —> J/ip + 7 . At first sight, this 
is surprising because the Lorentz structure of the helicity density matrix elements of XcJ 
production look more complicated 38, |39]. b However, after summing over the polarization 


J z of the XcJ meson, the helicity density matrix elements may indeed be expressed in the 
form of Eq. ([ 6 ]) as well. To ensure that our simplified method is correct, we compared with 
previous LO calculations 


38 


40] to find agreement. 


We generate the Feynman diagrams by using the FeynArts package 41] and calculate the 
amplitude squares with the help of the FeynCalc package 42]. As for direct J/ij) production, 
we reproduce the LO results for p xx ,(n) in Ref. 35] by setting q 2 = 0, while our results for 
the 0(v 2 ) corrections are new. As for J/-0 production via feed down from Xcj mesons, the 


LO helicity density matrix elements are found to agree with the results of Ref. 


after 


simplifying the latter, but are represented here for the first time in compact form, and the 
0{y 2 ) corrections are again new. We recover our results in Ref. 34] by summing over the 
helicity indices A and A'. I 11 the Appendix, we list all the new results for (n) and ( n ) 
in analytic form. 


III. PHENOMENOLOGICAL RESULTS 


We are now in a position to investigate the phenomenological significance of the 0{y 2 ) 
corrections to the polarization parameters in prompt J/'ijj photoproduction and hadropro- 
duction. I 11 our numerical analysis, we use m c = 1.5 GeV, a = 1/137.036, the LO 


In Eq. (22), the terms proportional to 


b We caution the reader of the following misprints in Ref. 

( O Xc0 ( 3 sf^))B(xcj —> J/ip+”f) with J = 1, 2 should be multiplied by the respective total-spin multiplicities 
(2J + 1). Three lines below Eq. (23), |Rp(0)| should be squared. The last factor in the first line of the 
expression for c-i in Eq. (A7) should read (7s 2 + st + It 2 ). 
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asymptotic scale parameter 


formula for aT s \p r ) with rij — 4 active quark flavors anc. 

Aqcd = 215 MeV 43], the CTEQ6L1 set for proton PDFs 43], the photon flux function 
given in Eq. (5) of Ref. 44] with ax = 2.5 GeV 2 6], and the choice pt r — P>f — y/Pr + 
for the renormalization and factorization scales. In contrast to the case of the unpolar¬ 
ized J/^} yield in Ref. |34|, the size of the 0(v 2 ) corrections to the J/%b polarization pa¬ 
rameters for a given cc Fock state n cannot be directly estimated from the SDC ratios 
R\\'(n ) = ( dF\\i(ri)/dx)/(dG\\i(n)/dx ), where x is some variable, because the dcw values 
enter as ratios in Eq. fl2]). For instance, if we had R 0 o(v) — Riiip) for some n, then the re¬ 
spective contribution to Xg would go unchanged upon the inclusion of the 0 {v 2 ) corrections. 
Therefore, we shall directly study the shifts induced in the polarization parameters by in¬ 
cluding the 0(v 2 ) corrections for each n separately. To standardize the numerical discussion, 
we quote the 0(v 2 ) corrections due to the 3 S'|^- 3 Df^ mixing relative to the corresponding LO 
results for n = 3 S ^. We exclude from our considerations the trivial J = 0 cases n = and 
n = 3 P^ ] . In want of fitted values of the LDMEs ( V H (n )), we estimate them with the help of 
the velocity scaling rule (V H (n)) / (O h (n)) = m 2 v 2 12] varying v 2 from 0 to 0.3 to account 
for uncertainties. Furthermore, we limit ourselves to the direct production of J/ip mesons 
and their production via the feed down from \cj mesons. In the latter case, we approxi¬ 
mately take into account the kinematic effect on the transverse momentum p t of the J/-0 
meson by setting pr = p^ cj Mj/^,/M XcJ , with Mj = 3.097 GeV, M Xcl = 3.511 GeV, and 


M Xc2 = 3.556 GeV 


45], which is justified since px M XcJ —Mj in typical experimental sit¬ 


uations. In the case of yp collisions, we only consider direct photoproduction for illustration 
because resolved photoproduction shares the partonic subprocesses with hadroproduction, 
which is studied separately. 

nn 

The J/ij} polarization was measured in a number of experiments 2-10]. As an illustration, 
we consider here three representative cases, namely HERA Run II [6], Tevatron Run II 
B, and the CMS experimental setup at the LHC 9], limiting ourselves to the hclicity 
frame. In HERA Run II 0. the polarization was measured in prompt photoproduetion 
at center-of-mass (c.m.) energy \fS = 319 GeV differential in px and inelasticity z = 
Pj/ii'Pp/P'y 'Ppi where p 7 , p p , and pj/^ are the four-momenta of the photon, proton, and J/%b 
meson, respectively, imposing in turn the acceptance cuts 0.3 < z < 0.9 and ]> 2 r > 1 GeV 2 
always in combination with the acceptance cut 60 GeV < W < 240 GeV on the yp c.m. 
energy W = (p 7 + p p ) 2 [6]. The 0(y 2 ) corrections to the polarization parameters A g, A^, 
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Pj-tGeV) Pr(GeV) 



(c) 


(d) 


FIG. 1: 0(u 2 )-corrected polarization parameters A g (solid lines), A^ (dotted lines), and \q$ (dashed 
lines) in direct J/tjj photoproduction at HERA Run II 6| through the cc Fock states (a) n = , 

(b) 3 sf, (c) 3 Pf\ and (d) 3 S'j 8 L 3 Dj 8 ] mixing as functions of px- The true results should lie inside 
the bands encompassed by the evaluations with v 2 = 0.3 (thick lines) and the LO results, with 
i; 2 = 0 (thin lines). 


and A e<j> in direct photoproduction at HERA Run II [6] through the cc Fock states n = 3 5'[ 1 ’ 8 ^, 
3 Pf \ and 3 5'[ 8 ' '’Df^ mixing are shown in Figs. |Tj and [2] as functions of px and z, respectively. 
The true results should lie inside the bands encompassed by the evaluations with v 2 = 0.3 
and the LO results, corresponding to v 2 = 0. We observe from Fig. |Tj that the 0(v 2 ) 
corrections are generally more significant in the small-pr range, while they tend to fade out 
for asymptotically large values of pr, with the exception of A^ for n = 3 5'j*' and of A g for 
3 S?- 3 D? mixing. Because of the steep falloff of the px distributions, the z distributions in 
Fig. [2] receive dominant contributions from the px region close to the cut at 1 GeV, toward 
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FIG. 2: Same as in Fig. [I] 
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(b) 
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(d) 

but as functions of z. 


the left ends of the panels in Fig. CD We observe from Fig. [2] that the 0{y 2 ) corrections to 
\e and Xg^ tend to be more sizable toward large values of z. The 0(v 2 ) corrections to A^ 
are relatively modest, except for the case of 3 ,S'{ 8 ^- 3 d[ 8 ^ mixing, where they are appreciable 
in the vicinity of z = 0.7. 

The px dependence of the polarization in prompt J/tjj hadroproduction was measured in 
Tevatron Run II at \fS = 1.96 TeV for \y\ < 0.6 3] and under CMS experimental conditions 
at yfS = 7 TeV for \y\ < 2.4 3. In both cases, we exclude from our considerations 
the small-pT’ range, Pt < 3 GeV, where the application of fixed-order perturbation theory 
is problematic. The 0[y 2 ) corrections to the polarization parameters Xg, A^, and A^ in 
prompt hadroproduction through the cc Fock states n = , 3 p[ 1 ’ 8 ^, and 3 S , j 8 ' :i D^ mixing 


are shown as functions of px for Tevatron Run II p] and the CMS setup 9] in Figs. [3] and QJ 
respectively. We observe from Fig. [3] that the 0{v 2 ) corrections to A g for n = 3 5f^, :i P\ l \ 
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are significant at small values of px, but quickly fade out toward large values of px, while 
they increase with px for n = 3 Pf ^ and :i S\ 8 ^ mixing. On the other hand, the 0(v 2 ) 
corrections to A^, are generally small and fade out as the value of px increases, except for 
the case of n = where the 0 (v 2 ) correction tends toward a constant negative value for 
increasing value of px, while the LO result tends to zero. Finally, the 0{y 2 ) corrections to 
A etj, are generally very small for all values of px■ The situation is very similar for the CMS 
setup [9], as may be seen by comparing Figs. [3] and [4] 


IV. SUMMARY 


In this work, we systematically calculated the relativistic corrections of relative order 
0(v 2 ) to the polarization observables Xg, A0, and A e<f> in prompt J/-0 photoproduction and 
hadropro duct ion, including both the direct mode and the feed down from XcJ an d mesons, 
and provided in analytic form the SDCs of all the relevant partonic subprocesses, which are 
listed in Eq. ([5]). This study is a natural extension of the one presented in Ref. 34], where the 
0{v 2 ) corrections to the unpolarized J/^ yields in photoproduction and hadroproduction 
were considered. All these 0(v 2 ) corrections are to be included along with the respective 
0(a s ) corrections 14-123] to complete the NLO treatments. We obtained the \ c j feed-down 


contributions in a compact form with the structure that is farni 


We found agreement with the LO results for direct production [33] and \ c j feed down 


iar from direct production 


38 


40]. 


Our results for the 0(v 2 ) corrections are new. Upon summation over the hclicity labels, we 
recover from them our Oiv 2 ) results for the unpolarized yield of prompt J/ip production 34], 
We numerically estimated the 0(v 2 ) corrections to the J/^ polarization parameters A g, A^, 
and Xg^ in direct photoproduction at HERA 6] and prompt hadroproduction at the Tevatron 
[3] and the LHC ijj] due to the relevant cc Fock states assuming that their LMDEs obey 


the velocity scaling rules 12]. We found that the shifts in Xg, X and A g^ due to the 0{v 2 ) 
corrections are significant in most cases, typically reaching ±0.2. Due to their nontrivial 
dependencies on p T and 2 and the characteristic differences between photoproduction and 
hadroproduction, their inclusion in global data analyses of J/ijj yield and polarization, in 
addition to the available 0 (a s ) corrections 14t23], may help to improve the quality of the 
determinations of the CO LDMEs and hopefully to remedy the J/-0 polarization crisis. 
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FIG. 3: 0(u 2 )-corrected polarization parameters Xg (solid lines), (dotted lines), and Xg$ (dashed 
lines) in prompt J/ijj hadroproduction at Tevatron Run II [3j through the cc Fock states (a) 
n = 3 SW, (b) 3 S^\ (c) 3 Pf\ (d) 3 5[ h L 3 Z?[^ mixing, (e) 3 p|^, and (f) 3 P^ as functions of pr- The 
true results should lie inside the bands encompassed by the evaluations with v 2 = 0.3 (thick lines) 
and the LO results, with v 2 = 0 (thin lines). 
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(e) (f) 

FIG. 4: Same as in Fig. 02 but for the CMS experimental setup at the LHC 
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Appendix A: Analytic results 


In this Appendix, we present in compact analytic form the coefficients A 2 £ (n) and B^ (n) 
dehned in Eq. (JH]) for the partonic subprocesses in Eq. ((5]). We refrain from presenting 
A^(n) for n = 3 S'[ 1 ’ 8 ^, 1 S'q 8 ^ 3 Pj 8 ^ because these expressions may be found in Appendix B 


of Ref. 


35j |. We list A l] (n) for n = 3 P\ X \ and B l] (n) for n = 3 Pj° J , 'PIO, and 3 Sf 


3chi8] 3 pi®] 3p[0 


3 cl®! 


3 .D® mixing. It is convenient to pull out common factors by writing A 2 ^ (n) = a* J (n).r//. and 
B 2 ^ (n) = b^(n)hk. The Mandelstam variables s, t, and u appearing below are dehned as 
s — (fci + fc 2 ) 2 , t — (k\ — P) 2 | q=0 , and u = (k 2 — P) 2 | q =o, and satisfy s + t + u = 4 m 2 . 


1. Photoproduction 


9 + 1 -A cc( 3 SP) + g: 


hi 


VnfsW) = _ 40967r 3 aag _ 

! 729 m c (4 m? c — s ) 3 (4m;? — t ) 3 (s + t ) 3 


(Ala) 


2048m 30 (3s 2 + 2 st + 3t 2 ) - 256m 8 (5s 3 - 12s 2 t + 5t 3 ) - 64m 6 (15s 4 + 80s 3 t 
90s 2 t 2 + 44st 3 + 15f 4 ) + 16m 4 (21s 5 + 103s 4 f + 182s 3 f 2 + 158s 2 t 3 + 67sf 4 + 21t 5 ) 
4m 2 (7s 6 + 42s 5 t + 99s 4 f 2 + 124s 3 t 3 + 87s¥ + 30sf 5 + 7t 6 ) 

3st(s + t) (s 2 + st + t 2 ) 2 (Alb) 


h 2 = 8 4096m 32 - 768m 30 (s - 3 1) + 64m 8 (5s 2 - 50st - 43f 2 ) + 16m 6 (8s 3 + 67 s 2 t 
+ 98 st 2 + 47f 3 ) + 4m 4 (2s 4 - 20s 3 t - 59 sH 2 - 58 st 3 - 25t 4 ) + m 2 t(- 6 s 4 - 4s 3 t 


2\2 


+ 15sV + 21st 3 + 12t 4 ) - t 2 (s 2 + st + t 2 ) 


(Ale) 


h 3 = 8 512m 8 (s 2 + t 2 ) + 16m 6 (13s 3 + 33s 2 t + 9sf 2 + 13f 3 ) - 4m 4 (13s 4 + 58 s 3 t 
+ 74 s 2 t 2 + 46 st 3 + 25f 4 ) + m 2 (s + t) (12s 4 + 29s 3 t + 48s 2 f 2 + 29 st 3 + 12f 4 ) 


2\2 


- (s + t ) 2 (s 2 + st + t 2 ) 


(Aid) 


h 4 = 8 128m 8 (5s 2 — t 2 ) + 16 m 6 (5s 3 + 19s 2 f + 14sf 2 + 8 f 3 ) + 4m 4 (s 4 — 22 s 3 t 
- 49 s 2 t 2 - 36 st 3 - 14t 4 ) + m 2 c t (7s 4 + 29s 3 t + 48s 2 t 2 + 34sf 3 + 12t 4 ) 


2\2 


— t(s + t) (s 2 + st + t 2 ) 


(Ale) 
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9 + 1 cc( 3 S) 8] ) + 3: 


^7(3^8]) = 1^97 ^[l] 
8 

The coefficients h *. are the same as those in Eq. (1 A1 [) . 
3 + 7 cc( 3 Pj 8] ) + g: 


V”( a P l f ] ) = 


1287r 3 aa^ 


15m 3 (s — 4m;?) 4 (£ — 4m;?) 4 ( s _|_ £)4 (_4 m 2 _|_ s _|_ £) 2 


(A2) 


(A3a) 


hi = (-4 ml + s + t) - 524288m 38 (s + t) (3s 2 + 14st + 3t 2 ) + 65536m 36 (3s - t)(3s 3 

- 3 sH + 25 st 2 - 9t 3 ) + 16384m 34 (39s 5 + 421s 4 t + 568s 3 t 2 + 908s 2 t 3 + 521st 4 + 39t 5 ) 

- 4096m 32 (126s 6 + 1064s 5 f + 2089s 4 t 2 + 3138s¥ + 2909sV + 1064st 5 + 126t 6 ) 

+ 1024m 30 (144s 7 + 1050s 6 t + 2473s 5 t 2 + 4177s¥ + 4857s¥ + 3273s¥ + 930st 6 
+ 144t 7 ) - 256m 3 (75s 8 + 349s 7 t + 816s 6 t 2 + 1732s 5 * 3 + 2564s 4 t 4 + 2432s 3 * 5 

+ 1376s 2 * 6 + 249s* 7 + 75t 8 ) + 64m 6 (s + t)(15s 8 - 66 s 7 1 - 314s 6 * 2 - 374s 5 * 3 

- 230s 4 * 4 - 254s 3 * 5 + 26s 2 * 6 - 106s* 7 + 15* 8 ) + 16m 4 s*(s + *)(39s 7 + 148s 6 * 

+ 219s 5 * 2 + 193s 4 * 3 + 213s 3 * 4 + 199s 2 * 5 + 88 st 6 + 39t 7 ) + 4m 2 s 2 * 2 (s + *)(13s 6 + 90s 5 * 

+ 232s 4 * 2 + 286 s 3 t 3 + 212sY + 70s* 5 + 13* 6 ) - lls¥(s + *) 2 (s 2 + st + t 2 ) 2 1 (A3b) 


h 2 = -16m 2 - 2097152m 20 (s + t) + 131072m 38 (31s 2 - 70s* - 5* 2 ) 

- 65536m 36 (l8s 3 - 191s 2 * - 180st 2 - 61* 3 ) - 8192m 34 (26s 4 + 1098s 3 * + 1287s 2 * 2 
+ 776 st 3 + 349t 4 ) + 4096m 32 (42s 5 + 775s 4 * + 1084s 3 * 2 + 556s 2 * 3 + 238s* 4 + 161* 5 ) 

- 512m 30 ( 79s 6 + 972s 5 * + 1074s 4 * 2 - 584s 3 * 3 - 1622s 2 * 4 - 1350s* 5 - 209* 6 ) 

+ 256m 8 (l6s 7 + 19s 6 * - 602s 5 * 2 - 1972s 4 * 3 - 2788s 3 * 4 - 2582s 2 * 5 - 1612st 6 - 365* 7 ) 

- 32/n 6 (2s 8 - 246s 7 * - 1554s 6 * 2 - 4296s 5 * 3 - 6520s 4 * 4 - 6780s 3 * 5 - 5509s 2 * 6 

- 3008s* 7 - 673* 8 ) - 16m 4 *(36s 8 + 205s 7 * + 664s 6 * 2 + 1260s 5 * 3 + 1731s 4 * 4 + 1755s 3 * 5 
+ 1412s 2 * 6 + 702s* 7 + 143* 8 ) + 2m 2 * 2 (-26s 8 - 54s 7 * + 5s 6 * 2 + 228s 5 * 3 + 648s 4 * 4 

+ 790s 3 * 5 + 669s 2 * 6 + 296s* 7 + 48* 8 ) - s* 3 (s + *)(6s 6 + 21s 5 * + 34s 4 * 2 + 42s 3 * 3 
+ 32s 2 * 4 + 19s* 5 + 6* 6 )1 (A3c) 
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h 3 = -16 m 2 c [262144m 46 (s + t) ( s 2 + t 2 ) - 8192m 44 (21s 4 + 218s 3 t + 114s 2 t 2 + 258st 3 
+ 21t 4 ) + 4096m 42 (9s 5 + 210s 4 t + 328s 3 t 2 + 648s 2 t 3 + 295st 4 + 14t 5 ) 

+ 2048m 40 (34s 6 + 72s 5 t + 50s 4 t 2 - 359s 3 t 3 - 490sV - 33st 5 + 14t 6 ) 

- 256m 3 (s + t) (218s 6 + 821s 5 t + 1234s 4 t 2 + 605s 3 t 3 - 31s 2 t 4 + 486st 5 + 143t 6 ) 

+ 32 m 6 (s + t)(497s 7 + 2427s 6 t + 4644s 5 t 2 + 5088s 4 t 3 + 3648s 3 t 4 + 2494sV 
+ 1697st 6 + 397t 7 ) - 16m 4 (s + t)(125s 8 + 745s 7 t + 1732s 6 t 2 + 2347s 5 t 3 + 2260s 4 t 4 
+ 1767s 3 t 5 + 1162s 2 t 6 + 575st 7 + 115t 8 ) + 2m 2 (s + t) 2 (48s 8 + 312s 7 t + 725s 6 t 2 
+ 970s 5 t 3 + 938s 4 t 4 + 770s 3 t 5 + 505s 2 t 6 + 252st 7 + 48t 8 ) - st(s + t) 3 (6s 6 + 21s 5 t 
+ 38sV + 48sV + 38sV + 21st 5 + 6t 6 )] (A3d) 


h 4 = -16m 2 [l048576m 48 (s - t)(s + t) - 98304m 46 (5s 3 - list 2 - 6 t 3 ) 

+ 8192m 44 (s 4 - 255s 3 t - 197s 2 t 2 - 93st 3 + 8 t 4 ) + 2048m 42 (90s 5 + 897s 4 t 
+ 1281s 3 t 2 + 645s 2 t 3 - 87st 4 - 122t 5 ) - 512m 40 (226s 6 + 1467s 5 t + 2568s 4 t 2 
+ 1908s 3 t 3 + 262s 2 t 4 - 895st 5 - 3721 6 ) + 128m 8 (223s 7 + 1197s 6 t + 2055sV 
+ 1372s 4 t 3 - 390s 3 t 4 - 1905s 2 t 5 - 2010st 6 - 614t 7 ) - 32 m 6 (99s 8 + 404s 7 t + 282s 6 t 2 

- 989s 5 t 3 - 2828s 4 t 4 - 3803s s t 5 - 3693s 2 t 6 - 2248st 7 - 556t 8 ) + 8m 4 (18s 9 + 16s 8 t 

- 385s 7 t 2 - 1503s 6 t 3 - 2951s 5 t 4 - 3933s 4 t 5 - 3790s 3 t 6 - 2794s 2 t 7 - 1296st 8 - 2581 9 ) 
+ 2 m 2 t(s + t)(16s 8 + 145s 7 t + 409s 6 t 2 + 681s 5 t 3 + 877s 4 t 4 + 836s 3 t 5 + 602s 2 t 6 

+ 274st 7 + 48t 8 ) - st 2 (s + t ) 2 (7s 6 + 23s 5 t + 39s 4 t 2 + 47s 3 t 3 + 36s 2 t 4 + 20st 5 + 6 t 6 ) 

(A3e) 

9 + 7 ^ cc( 3 Sf\ 3 Df ] ) + g: 


b 9 \ 3 SP, 3 Df l ) = 


17r 3 aa 2 


81m c (s — 4 m;?) 4 (t — 4 m;?) (s + t) 4 


(A4a) 


hi — — (4m 2 — s) (4m 2 — t) (s + t) 2048m 40 (3s 2 + 5st + 3t 2 ) — 256m 8 (5s 3 + 24s 2 t 
+ 24st 2 + 5t 3 ) - 64m 6 (15s 4 + 2s 3 t - 18s 2 t 2 + 2st 3 + 15t 4 ) + 16m 4 (21s 5 + 31s 4 t 
+ 20s 3 t 2 + 20s 2 t 3 + 31st 4 + 21t 5 ) - 4m 2 (7s 6 + 12s 5 t + 9s 4 t 2 - 2s 3 t 3 + 9s 2 t 4 + 12st 5 
+ 7t 6 ) — 3st(s + t) (s 2 + st + t 2 ) 2 (A4b) 


16 



h 2 = -4 (4 ml - t) [81920m 44 (s + t) - 2048m 42 (25s 2 + 52st + 15t 2 ) 

+ 512m 40 (26s 3 + 75s 2 t + 36st 2 - 13t 3 ) - 128m 3 (3s 4 + 18s 3 t - 24s 2 t 2 - 92st 3 
- Alt 4 ) - 32m 6 (l2s 5 + 39s 4 t + 100s 3 t 2 + 184s 2 t 3 + 166st 4 + 43t 5 ) + 8m 4 (4s 6 
+ 30s 5 t + 81s 4 t 2 + 140s 3 t 3 + 161s 2 t 4 + 106st 5 + 30t 6 ) - 2m 2 t(6s 6 + 12s 5 t + 29s 4 t 2 
+ 38s 3 t 3 + 41s 2 t 4 + 24st 5 + 10 t 6 ) — st 2 (s + t ) (s 2 + st + t 2 ) 2 (A4c) 


, 3 = —4(s + t) [4096m 42 (s 2 + t 2 ) + 512m 40 (29s 3 + 61s 2 t + 61st 2 + 29t 3 ) 

- 256m 3 (s + t ) 2 (37s 2 + 43st + 371 2 ) + 32m 6 (73s 5 + 299s 4 t + 536s 3 t 2 + 536s 2 t 3 
+ 299st 4 + 731 5 ) - 16m 4 (20s 6 + 80s 5 t + 167s 4 t 2 + 212s 3 t 3 + 167s 2 t 4 + 80st 5 + 20t 6 ) 

+ 2m 2 (10s 7 + 32s 6 t + 59s 5 t 2 + 77s 4 t 3 + 77s 3 t 4 + 59s 2 t 5 + 32st 6 + 10t 7 ) + st(s 3 
+ 2 s 2 t + 2 st 2 + t 3 ) 2 ] (A4d) 

h 4 = -4 2048m 42 (l3s 3 + 7s 2 t - 5st 2 - lit 3 ) - 512m 40 (14s 4 + 5s 3 t - 61s 2 t 2 - 95st 3 
- 43t 4 ) - 128m 3 (s 5 + 43s 4 t + 214s 3 t 2 + 340s 2 t 3 + 265st 4 + 731 5 ) + 32m 6 (2s 6 
+ 71s 5 t + 305s 4 t 2 + 536s¥ + 530s 2 t 4 + 301st 5 + 71t 6 ) - 8m 4 t(40s 6 + 167s 5 t 
+ 345s 4 t 2 + 413s 3 t 3 + 327s 2 t 4 + 160st 5 + 40t 6 ) + 2m 2 t(10s 7 + 30s 6 t + 57 s 5 t 2 
+ 77s 4 t 3 + 79s 3 t 4 + 61s 2 t 5 + 32st 6 + 10 t 7 ) + st 2 (s 3 + 2 sH + 2 st 2 + t 3 ) 2 1 (A4e) 


q{q) + 7 ^ cc( 3 Af ] ) + q(q): 


„ roi Air^aalel 

59(9)7 pcd 8 ! \ __ s g 

4 27m®st (Ami ~ s ) 


(A5a) 


hi = —640m 6 + 160m 4 (2s + t) — 4m 2 (27s 2 + lOst + 51 2 ) + 11s (s 2 + t 2 ) (A5b) 


h 2 = —16m 2 (44m 2 — 5s) 


(A5c) 


h 3 = —32m 2 (44m 2 — 5s) 


(A5d) 


q{q) + 7 -t cc( 3 Pj 3] ) + g(g): 


59(9)7 ( 3 pj 3 ]) = 


h A = —16m 2 (44m 2 — 5s) 


2567r 3 aa 2 

135m® (Am 2 c — s) (s + t) 4 (— Am 2 c + s + t) z 


(A5e) 


(A6a) 


17 



hi = (s + t) (s + t — 4m 2 ) 1280m® (s + 5t) — 128?n^t(6s + 5t) + 4m 2 (21s 
+ 49s 2 t + st 2 + 21t 3 ) — 11s (s 3 + s 2 t + st 2 + t 3 ) 


(A6b) 


hi = -32 ml 


1024m 3 - 128m® (s - 14t) + 16m* (21s 2 + 14st - 29 1 2 ) - 4m 2 (4s 3 


+ 13s 2 t + 54st 2 + 3t 3 ) + s (s 3 + s 2 t + 3st 2 + 3t 3 ) 


(A6c) 


h 3 = —64m 2 (s + t) 80m^(3s + 7 1) — 8m 2 (3s 2 + list + 8t 2 ) + s(s + t ) 2 (A6d) 


hi = —64m 2 64m®(s — t) + 8m* (19s 2 + 72 st + 191 2 ) — 2m 2 (s 3 + 27 s 2 t + 43st 2 
+ 17t 3 ) + st(s + t) 2 (A6e) 


q(q) + 7 -> cc^f 1 , 3 H[ 81 ) + g(g): 


69 ( 9 ) 7 ( 35 { 8] j3 £>[ 8 ]) = 


8 \ / |7r 3 aa^eg 

9 m^st 


hi = 32 m* — 8 m 2 (s + i) + s 2 + t 2 


/12 = 16 m 2 
h 3 = 32 m 2 
h 4 = 16m 2 

2. Hadroproduction 

3 + 3 -+ cc( 3 5l 1] ) + 3: 

&"( 3 4 1] ) = |^fe 97 ( 3 Al 1] ). 

The coefficients hk are the same as those in Eq. (I A II) . 


(A7a) 

(A7b) 

(A7c) 

(A7d) 

(A7e) 


(AS) 
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9 + 9 -»• cc( 3 Pj 1] ) + g: 



I67r 3 a 3 

9m 3 (s — Ami) 4 (t — Am 2 c ) 4 (s + t) 4 


(A9a) 


gi = (s — Ami) (4' m c — t) (s + t) [16m 4 (s 2 + st + t 2 ) — 4m 2 (2s 3 + 3s 2 t + 3st 2 
+ 2t 3 ) + (s 2 + st + t 2 y ] 16m®(s + t) — Ami (5-s 2 + 6 st + 5t 2 ) + m 2 (8s 3 + 13s 2 t 
+ 13st 2 + 8t 3 ) — (s 2 + st + t 2 y (A9b) 


<72 = —4m 2 (t —4m 2 ) 2 512m 40 (s + t) 2 — 768m 3 (s + t) 3 + 32m 6 (13s 4 + 56s 3 t 
+ 86s 2 t 2 + 58st 3 + 15t 4 ) - 16m 4 (6s 5 + 40s 4 t + 82s 3 t 2 + 85s 2 t 3 + 45st 4 + 10t 5 ) 

+ 2 m 2 c (4s 6 + 54s 5 t + 137s 4 t 2 + 192s 3 t 3 + 153sY + 70st 5 + 14t 6 ) 

— t (s 2 + st + t 2 ) 2 (7s 2 + 7st + 2t 2 ) (A9c) 


r/ 3 = 4m 2 (s + t) 2 512m 40 (s 2 + t 2 ) — 256m 3 (3s 3 + 2s 2 t + 2st 2 + 3t 3 ) 

+ 32m 6 (l5s 4 + 18s 3 t + 10s 2 t 2 + 18st 3 + 15t 4 ) — 16m 4 (10s 5 + 15s 4 t + 9s 3 t 2 + 9 s 2 t 3 
+ 15st 4 + 10t 5 ) + 2m 2 (14s 6 + 22s 5 t + 17s 4 t 2 + 8s 3 t 3 + 17s 2 t 4 + 22st 5 + 14t 6 ) 

— (s 2 + st + t 2 ) 2 (2s 3 — s 2 t — st 2 + 2t 3 ) (A9d) 


r/ 4 = 4m 2 (s — 4m 2 ) 512m 40 (s — t)(s + t) 2 — 128m® (7s 4 + 14s 3 t + 2s 2 t 2 — 12st 3 
- 7t 4 ) + 32m 6 (18s 5 + 53s 4 t + 46s 3 t 2 - 18s 2 t 3 - 44st 4 - 19t 5 ) - 8m 4 (22s 6 + 86s 5 t 
+ 119s 4 t 2 + 46s 3 t 3 - 60s 2 t 4 - 72st 5 - 251 6 ) + 2m 2 (14s 7 + 64s 6 t + 114s 5 t 2 + 95s 4 t 3 
+ 3s 3 t 4 — 63s 2 t 5 — 55st 6 — 16t 7 ) — (s 2 + st + t 2 ) 2 (2s 4 + 5s 3 t + 3s 2 t 2 — 4st 3 — 2t 4 ) 

(A9e) 

ff 9 ( 3 p[ 1] ) = ----g- (A9f) 

45m 3 {Am 2 c — s) 5 {Am 2 c — t) 5 (s + t) 5 
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hi = (4m 2 — s) (4m 2 — t) (s + t ) 65536m 36 (s 3 + s 2 t + st 2 + t 3 ) — 4096m 34 (7s 4 

- 73s 3 t - 60s¥ - 53st 3 + 7t 4 ) - 1024m 32 (132s 5 + 853s 4 t + 1235s¥ + 1115s¥ 

+ 573st 4 + 92t 5 ) + 256m 30 (686s 6 + 3499s 5 t + 6337s 4 t 2 + 6952s¥ + 5057sV 

+ 2179st 5 + 446t 6 ) - 64m 3 (1420s 7 + 7165s 6 t + 15279s¥ + 20080s¥ + 18020s¥ 

+ 11059s¥ + 4205st 6 + 860t 7 ) + 16m 6 (1485s 8 + 7888 s 7 t + 19416s 6 t 2 + 30007s 5 t 3 
+ 32040s 4 t 4 + 24567s 3 t 5 + 13136s 2 t 6 + 4428st 7 + 845t 8 ) - 4m 4 (780s 9 + 4522s 8 t 
+ 12768s 7 t 2 + 22889s 6 t 3 + 28781s 5 t 4 + 26401s 4 t 5 + 17649s¥ + 8348sV + 2482st 8 
+ 420t 9 ) + m 2 (s 2 + st + t 2 ) 2 (164s 6 + 792s 5 t + 1609s 4 t 2 + 1762s¥ + 1289s 2 t 4 
+ 472st 5 + 84t 6 ) — 11 st(s + t) (s 2 + st + t 2 ) 4 (A9g) 

h 2 = 4m 2 (t - 4m 2 ) 2 [8192m 34 (s + t) 2 (19s + 39 1) - 2048m 32 (161s 4 + 696s 3 t 
+ 1206s 2 t 2 + 928 st 3 + 257t 4 ) + 512 m 30 (473s 5 + 2320s 4 t + 5423s¥ + 6065s 2 t 3 
+ 3340st 4 + 7231 5 ) - 128m 8 (623s 6 + 3752 s b t + 11313s 4 t 2 + 16704s¥ + 13939s 2 t 4 
+ 6152st 5 + 1149t 6 ) + 32m 6 (376s 7 + 3161s 6 t + 11876s 5 t 2 + 22019s¥ + 24449s¥ 

+ 16498s¥ + 6351st 6 + 1086t 7 ) - 8m 4 (84s 8 + 1328s 7 t + 6071s¥ + 13796s 5 t 3 
+ 18908s 4 t 4 + 17172 s¥ + 10043s 2 t 6 + 3536st 7 + 566t 8 ) + 2m 2 t(218s 8 + 1208s 7 t 
+ 3361s 6 t 2 + 5431 s¥ + 6167s 4 t 4 + 4797s 3 t 5 + 2558s¥ + 832st 7 + 124t 8 ) 

- st 2 (s 2 + st + t 2 ) 2 (7s 3 + 14s 2 t + 9st 2 + 2t 3 ) (A9h) 

h 3 = —4m 2 (s + t) 2 262144m 36 (s 2 + t 2 ) — 8192m 34 (9s 3 — 19s 2 t + 21st 2 + 9t 3 ) 

- 2048m 32 (157s 4 + 400s 3 t + 218s 2 t 2 + 120st 3 + 137t 4 ) + 512m 30 (683s 5 + 1914s 4 t 
+ 1569s 3 t 2 + 849s 2 t 3 + 994st 4 + 563t 5 ) - 128m 8 (1321s 6 + 4100s 5 t + 4181s 4 t 2 

+ 2372s¥ + 2001s 2 t 4 + 2340st 5 + 1021t 6 ) + 32m 6 (1382s 7 + 4555s 6 t + 5490s 5 t 2 
+ 3491s¥ + 2011s 3 t 4 + 2610s¥ + 2595st 6 + 1022t 7 ) - 8m 4 (750s 8 + 2584s 7 t 
+ 3563s¥ + 2580s 5 t 3 + 1242s¥ + 1100s 3 t 5 + 1703s 2 t 6 + 1424st 7 + 550t 8 ) 

+ 2 m 2 (164s 9 + 600s 8 t + 926s 7 t 2 + 753s 6 t 3 + 383s 5 t 4 + 183s¥ + 313s s t 6 + 446sV 
+ 320st 8 + 124t 9 ) - st (2s 2 - 3st + 21 2 ) (s 3 + 2s 2 t + 2st 2 + t 3 ) 2 (A9i) 
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hi = 4 m 2 (4 m 2 — s) 8192-m 14 (13s 4 + 76s 3 t + 42 s 2 t 2 — 52 st 3 — 311 4 ) 

- 4096m 12 (92s 5 + 435s 4 t + 486s¥ + 13s 2 t 3 - 322st 4 - 128t 5 ) + 512m 40 (729s 6 

+ 3578 s 5 t + 6026s 4 t 2 + 3382 s 3 t 3 - 1629s¥ - 2648st 5 - 846t 6 ) - 256m®(664s 7 

+ 3613s 6 t + 7719s 5 t 2 + 7378 sH 3 + 1751s¥ - 2993 s 2 t 5 - 2634 st 6 - 722 1 7 ) 

+ 32m 6 (1324s 8 + 7870 s 7 t + 19424s 6 t 2 + 24728s 5 t 3 + 14981s 4 t 4 - 1074s 3 t 5 

- 8562 s 2 t 6 - 5580 st 7 - 1351t 8 ) - 16m 4 (359s 9 + 2242s 8 t + 6074s 7 t 2 + 9229s 6 t 3 
+ 7941s 5 t 4 + 2897s 4 t 5 - 1764s 3 * 6 - 2790s 2 * 7 - 1522st 8 - 330t 9 ) + 2m 2 (164s 10 

+ 1060s 9 t + 3088s 8 * 2 + 5238 s 7 t 3 + 5563s 6 * 4 + 3338s 5 * 5 + 368s 4 * 6 - 1476s 3 * 7 - 1411s 2 * 8 

- 680s* 9 — 1327 10 ) — st (s 2 + st + t 2 ) 2 (2s 5 + 7s 4 * + 8s 3 * 2 — s 2 t 3 — 6st 4 — 2t 5 ) (A9j) 


g + g ->■ cc( 3 pj 1] ) + g: 


a fl , 9('3p[l]\ __ a s _ 

135 m^st (s — 4m^) 4 (* — 4-m^) 4 (s + *) 4 (—4m^ + s + t) 


(AlOa) 


gi = -18874368m 20 (s + t) 4 + 9437184m 18 (s + t) 3 (3s 2 + 5st + 3t 2 ) 

- 16384m 16 (s + t) 2 (1224s 4 + 3897s 3 t + 5297 s 2 t 2 + 3897 st 3 + 1224* 4 ) 

+ 4096-m 14 (2160s 7 + 11913s 6 * + 29647s 5 * 2 + 44652s 4 * 3 + 44652s 3 * 4 + 29647s 2 * 5 
+ 11913s* 6 + 2160* 7 ) - 1024m 12 (2592s 8 + 14994s 7 * + 39656s 6 * 2 + 65365s 5 * 3 
+ 76238s 4 * 4 + 65365s 3 * 5 + 39656s 2 * 6 + 14994s* 7 + 2592* 8 ) + 256 m 10 (2160s 9 
+ 13374s 8 * + 37078s 7 * 2 + 63603s 6 * 3 + 79713s 5 * 4 + 79713s 4 * 5 + 63603s 3 * 6 + 37078s 2 * 7 
+ 13374s* 8 + 2160* 9 ) - 64m® (1224s 10 + 8445s 9 * + 24967s 8 * 2 + 43880s 7 * 3 + 55214s 6 * 4 
+ 58192s 5 * 5 + 55214s 4 * 6 + 43880s 3 * 7 + 24967s 2 * 8 + 8445s* 9 + 1224* 10 ) + 16m 6 (432s n 
+ 3501s 10 * + 11595s 9 * 2 + 21568s 8 * 3 + 26848s 7 * 4 + 27020s 6 * 5 + 27020s 5 * 6 + 26848s 4 * 7 
+ 21568s 3 * 8 + 11595s 2 * 9 + 3501st 10 + 432* 11 ) - 4-m 4 (72s 12 + 792 s u t + 3218s 10 * 2 
+ 6815s 9 * 3 + 8979s 8 * 4 + 8574s 7 * 5 + 7840s 6 * 6 + 8574s 5 * 7 + 8979s 4 * 8 + 6815s 3 * 9 
+ 3218s 2 * 10 + 792 st 11 + 72* 12 ) + m 2 c st (s 2 + st + t 2 ) 2 (60s 7 + 264s 6 * + 341s 5 * 2 + 71s 4 * 3 
+ 71s 3 * 4 + 341s 2 * 5 + 264st 6 + 60* 7 ) — s 2 t 2 (s + t) 2 (s 2 + st + t 2 ) 4 (AlOb) 
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g 2 = -12 m 2 (4 m 2 - t) [8192 m 44 (s + t ) 2 (12s 2 + 47 st + 24t 2 ) - 2048m 42 (56s 5 

+ 345s 4 t + 722 sh 2 + 717s 2 t 3 + 352st 4 + 681 5 ) + 512m 40 (112s 6 + 741s 5 t + 1745sV 
+ 2129s 3 t 3 + 1501sV + 552 st 5 + 72t 6 ) - 128m 8 (128s 7 + 919s 6 t + 2368s 5 t 2 
+ 3306s 4 t 3 + 2900s 3 t 4 + 1591s 2 t 5 + 480st 6 + 48t 7 ) + 32m® (92s 8 + 704s 7 t + 1933s 6 t 2 
+ 2902s 5 t 3 + 2844sV + 1945s¥ + 905s 2 t 6 + 259st 7 + 40t 8 ) - 8m 4 (40s 9 + 324s 8 t 
+ 952s 7 t 2 + 1461s 6 t 3 + 1416s 5 t 4 + 983s 4 t 5 + 558s 3 t 6 + 298 s 2 t 7 + 132st 8 + 36t 9 ) 

+ 2 m 2 (8s 10 + 80s 9 t + 274s s t 2 + 432sV + 355s 6 t 4 + 133s 5 t 5 + 5s 4 t 6 + 71sV + 132sY 
+ 94st 9 + 24t 10 ) - t (s 2 + st + t 2 ) 2 (4s 6 + 12s 5 t + sH 2 - 18s 3 t 3 - s 2 t 4 + 10st 5 + 4t 6 ) 

(AlOc) 

r /3 = — 12m 2 (s + t) [l572864m 46 (s + t) 3 - 294912m 44 (s + t ) 2 (7s 2 + lOst + 7 1 2 ) 

+ 2048m 42 (588s 5 + 2257s 4 t + 3917s 3 t 2 + 3917s 2 t 3 + 2257 st 4 + 5881 5 ) 

- 512m 40 (824s 6 + 3435s 5 t + 6696s 4 t 2 + 8202s 3 t 3 + 6696s 2 t 4 + 3435st 5 + 824t 6 ) 

+ 128 m 8 (784s 7 + 3461s 6 t + 7141s 5 t 2 + 9650s 4 t 3 + 9650sY + 7141s 2 t 5 + 3461st 6 
+ 784t 7 ) - 32m 6 (524s 8 + 2459s 7 t + 5234s 6 t 2 + 7181s 5 t 3 + 7796s 4 t 4 + 7181s 3 t 5 
+ 5234sV + 2459st 7 + 524t 8 ) + 8m 4 (236s 9 + 1218s 8 t + 2722 s 7 t 2 + 3705s 6 t 3 
+ 3905sV + 3905s 4 t 5 + 3705s 3 * 6 + 2722sY + 1218st 8 + 236t 9 ) - 2m 2 (s + t) 2 (64s 8 
+ 246s 7 t + 356s 6 * 2 + 337s 5 t 3 + 274s 4 t 4 + 337s s t 5 + 356s 2 t 6 + 246st 7 + 64t 8 ) 

+ (s 2 + st + t 2 y (4s 7 + 18s 6 t + 23 s 5 t 2 + 3s 4 t 3 + 3s 3 t 4 + 23 s 2 t 5 + 18st 6 + 4t 7 ) (AlOd) 
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g A = -12 m 2 [l96608m 16 (s + *) 3 (3s + 5 1) - 73728m 14 (s + *) 2 (12s 3 + 33s 2 * + 35s* 2 
+ 16* 3 ) + 2048m 32 (280s 6 + 1369s 5 * + 3023s 4 * 2 + 3917s 3 * 3 + 3151s 2 * 4 + 1476s* 5 

+ 308* 6 ) - 512m 10 (416s 7 + 2129s 6 * + 5074s 5 * 2 + 7460s 4 * 3 + 7438s 3 * 4 + 5057s 2 * 5 

+ 2130s* 6 + 408* 7 ) + 128 m 3 (394s 8 + 2169s 7 * + 5479s 6 * 2 + 8608s 5 * 3 + 9650s 4 * 4 
+ 8183s 3 * 5 + 5123s 2 * 6 + 2076s* 7 + 390* 8 ) - 32m 6 (244s 9 + 1504s 8 * + 4073s 7 * 2 
+ 6653s 6 * 3 + 7734s 5 * 4 + 7243s 4 * 5 + 5762s 3 * 6 + 3620s 2 * 7 + 1479s* 8 + 280* 9 ) 

+ 8m 4 (92s 10 + 674s 9 * + 2054s 8 * 2 + 3591s 7 * 3 + 4219s 6 * 4 + 3905s 5 * 5 + 3391s 4 * 6 
+ 2836s 3 * 7 + 1886s 2 * 8 + 780s* 9 + 144* 10 ) - 2m 2 (16s 11 + 162s 10 * + 606s 9 * 2 + 1224s 8 * 3 
+ 1555s 7 * 4 + 1400s 6 * 5 + 1126s 5 * 6 + 1044s 4 * 7 + 983s 3 * 8 + 680s 2 * 9 + 276s* 10 + 48* 11 ) 

+ * (s 2 + st + * 2 ) 2 (6s 7 + 23s 6 * + 28s 5 * 2 + 3s 4 * 3 — 2s 3 * 4 + 18s 2 * 5 + 16s* 6 + 4*') (AlOe) 

^ 2 ^ 675m^s 2 * 2 (4m;? — s) 5 (4m;? — *) 5 (s + *) 5 (—4m;? + s + *) 2 ^ ^ 
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h ) = -St (-Am 2 c + s + t) [301989888(s + t) 3 (l5s 2 + 38ts + 15t 2 )m 24 

- 8388608(s + t) 2 (945s 4 + 3943ts 3 + 5985 t 2 s 2 + 3763t s s + 945t 4 ) m 22 

+ 262144(27624s 7 + 179241ts 6 + 503190tV + 788233 t 3 s 4 + 755573t 4 s 3 + 450370tV 
+ 156201t 6 s + 24744t 7 )m 20 - 65536(65928s 8 + 438138ts 7 + 1296615t 2 s 6 
+ 2241746t 3 s 5 + 2524170tV + 1935826t 5 s 3 + 1003615t 6 s 2 + 325218t 7 s + 50088t 8 )m 48 
+ 16384(108144s 9 + 742491ts 8 + 2285700 t 2 s 7 + 4194310t 3 s 6 + 5180803tV 
+ 4615783t 5 s 4 + 3046670t e s 3 + 1457000tV + 456891t 8 s + 69264t 9 )m 46 

- 4096(124080s 10 + 900416ts 9 + 2906497t 2 s 8 + 5577856tV + 7244723t 4 s 6 

+ 6963392t 5 s 5 + 5278403t 6 s 4 + 3209416t 7 s 3 + 1491977t 8 s 2 + 460856t 9 s + 67920t 10 )m 44 
+ 1024(99240s 41 + 779395ts 10 + 2690450t 2 s 9 + 5443508t 3 s 8 + 7327604tV 
+ 7226103t 5 s 6 + 5774383t 6 s 5 + 4032804tV + 2433368t s s 3 + 1137870t 9 s 2 + 343315t 10 s 
+ 47400t n )m 42 - 256(53064s 12 + 463250ts n + 1757329tV° + 3842678t 3 s 9 
+ 5442560t 4 s 8 + 5431020t 5 s 7 + 4269078t 6 s 6 + 3095980t 7 s 5 + 2225920t s s 4 + 1409878t 9 s 3 
+ 655689t 10 s 2 + 185450t n s + 22824t 12 )m 40 + 64(16968s 13 + 173309ts 12 + 751960t 2 s n 
+ 1842678t 3 s 10 + 2850960t 4 s 9 + 2968913t 5 s 8 + 2256004t 6 s 7 + 1525084t 7 s 6 
+ 1187673t 8 s 5 + 968900t 9 s 4 + 624838t 10 s 3 + 270020t n s 2 + 66749t 12 s + 6888t 13 )m 8 

- 16(2424s 14 + 34060ts 13 + 183551t 2 s 12 + 532924t 3 s n + 949514tV° + 1102104t 5 s 9 
+ 872595t 6 s 8 + 540792t 7 s 7 + 397555t s s 6 + 391504t 9 s 5 + 330714t 10 s 4 + 190964t n s 3 
+ 68951t 12 s 2 + 13300t ls s + 984t 14 )m 6 + 4st(2116s 13 + 19052ts 12 + 74123t 2 s n 

+ 165228t 3 s 10 + 235243t 4 s 9 + 227226t 5 s 8 + 164836tV + 119456t 7 s 6 + 113946t s s 5 
+ 107443t 9 s 4 + 73328t 10 s 3 + 32643t n s 2 + 8372 t 12 s + 916t 13 )m 4 - s 2 t 2 (s 3 + 2 ts 2 
+ 2 t 2 s + t 3 ) 2 (360s 6 + 1356ts 5 + 1421t 2 s 4 + 494t 3 s 3 + llOltV + 1036t 5 s + 280t 6 ) m\ 

+ 7s 3 t 3 (s + t) 3 (s 2 + ts + t 2 ) 4 (AlOg) 
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h 2 = -12 m 2 c (4 ml - t) [l6777216st(s + tf (8s 2 - 9 ts + 2 1 2 ) mf + 524288(32s 7 

- 636 ts 6 - 2381t 2 s 5 - 4579t 3 s 4 - 5103t 4 s 3 - 2877t 5 s 2 - 608t 6 s + 32 t 7 )mf 

- 131072(224s 8 - 2256ts 7 - 13229t 2 s 6 - 34841t 3 s 5 - 50353t 4 s 4 - 39495tV - 16038t 6 s 2 

- 2292 t 7 s + 256t 8 )m 48 + 32768(672s 9 - 4044ts 8 - 34626t 2 s 7 - 111888t 3 s 6 

- 195319t 4 s 5 - 196987tV - 116723t 6 s 3 - 35861t 7 s 2 - 2592t 8 s + 896t 9 )m 46 

- 8192(1120s 10 - 4264ts 9 - 54962t 2 s 8 - 207478tV - 416241t 4 s 6 - 501995tV 

- 383465t 6 s 4 - 180583tV - 41448t 8 s 2 + 1468t 9 s + 1792t 10 )m 44 + 2048(1120s n 

- 3164ts 10 - 59449t 2 s 9 - 249681t 3 s 8 - 549542tV - 743186fV - 667344tV 

- 410272tV - 160085t 8 s 3 - 25309tV + 6312t 10 s + 2240£ n )m 42 - 512(672s 12 

- 2176ts n - 46457t 2 s 10 - 206373t s s 9 - 475554tV - 672814t 5 s 7 - 639440t 6 s 6 
_ 444400tV - 235223t 8 s 4 - 84109tV - 8818t 10 s 2 + 6012t n s + 1792t 12 )mJ° 

+ 128(224s 13 - 1396ts 12 - 26208 t 2 s u - 119918t 3 s 10 - 279940tV - 383196tV 

- 319618t 6 s 7 - 174546t 7 s 6 - 88697t 8 s 5 - 63753t 9 s 4 - 36485t 10 s 3 - 7707t n s 2 

+ 1808t 12 s + 896t 13 )m 8 - 32(32s 14 - 584ts 13 - 9872tV 2 - 47988tV 1 - 115308tV° 

- 146248t 5 s 9 - 68062t e s 8 + 62478tV + 108433t s s 6 + 47003t 9 s 5 - 20127t 10 s 4 

- 29125t n s 3 - 11012t 12 s 2 - 972t 13 s + 256t 14 )m 6 - 8t(104s 14 + 2056ts 13 + 11846tV 2 
+ 32068t 3 s n + 41516t 4 s 10 + 3659tV - 73197t 6 s 8 - 117854tV - 85086t s s 6 

- 15497tV + 23451t 10 s 4 + 20944t n s 3 + 7342t 12 s 2 + 952t 13 s - 32t 14 )m 4 

+ 2st 2 (144s 13 + 1304ts 12 + 4612t 2 s n + 7188t 3 s 10 + 183t 4 s 9 - 19901t 5 s 8 - 38638 t 6 s 7 

- 36862t 7 s 6 - 15497t 8 s 5 + 4991tV + 10688t 10 s 3 + 6400t n s 2 + 1860t 12 s + 232 t 13 )m 2 c 

- sH 3 (s 3 + 2 ts 2 + 2 t 2 s + t 3 ) 2 (4s 6 + 12ts 5 - 83t 2 s 4 - 186t 3 s 3 - 21t 4 s 2 + 74t 5 s + 20t 6 ) 

(AlOh) 
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h 3 = 12 m 2 c (s + t ) [l207959552st(s + t) 2 (2s 2 + 5 ts + 2t 2 ) mf - 2097152st(2003s 5 
+ 101791s 4 + 20112t 2 s 3 + 19932lV + 100891 4 s + 2093l 5 )m 20 - 131072(32s 8 

- 26276Is 7 - 146167l 2 s 6 - 337414l s s 5 - 428510lV - 324974l 5 s 3 - 140887l 6 s 2 

- 274761 7 s + 32t 8 )m 1 c 8 + 32768(224s 9 - 545881s 8 - 323187lV - 8236031 s s 6 

- 1219834l 4 s 5 - 1174194l 5 s 4 - 7478631 6 s 3 - 293687l 7 s 2 - 542281 8 s + 224l 9 )m 46 

- 8192(672s 10 - 801081s 9 - 497784l 2 s 8 - 1348997t 3 s 7 - 21856061 4 s 6 - 2420290fV 

- 19401661 6 s 4 - 1111197l 7 s 3 - 411864l s s 2 - 72268l 9 s + 672 t 10 }m l c A + 4096(560s n 

- 430421s 10 - 2819921 2 s 9 - 8004241 3 s 8 - 13619651 4 s 7 - 16142731V - 1466973lV 

- 10575251 7 s 4 - 577964l 8 s 3 - 207612lV - 344621 10 s + 5601 n )m 12 - 512(1120s 12 

- 676521s 11 - 4742511 2 s 10 - 14176881 3 s 9 - 25022611 4 s 8 - 30326981 5 s 7 - 28613161 6 s 6 

- 2322578l 7 s 5 - 163086UV - 889688l 9 s 3 - 31149U 10 s 2 - 482121 n s + 11201 12 )m 10 
+ 128(672s 13 - 380361s 12 - 288627l 2 s n - 922087lV° - 17080031 4 s 9 - 20998991V 

- 19336121 6 s 7 - 15745721 7 s 6 - 12757991V - 9400631V - 517827l 10 s 3 - 173847l n s 2 

- 246361 12 s + 6721 13 )m 8 - 32(224s 14 - 146281s 13 - 1202421V 2 - 4152691 3 s n 

- 821945lV° - 10454191 5 s 9 - 925223iV - 670140lV - 542503tV - 510619tV 

- 4065851 10 s 4 - 2201491 n s 3 - 690421 12 s 2 - 88681 13 s + 224l 14 )m 6 + 8(32s 15 

- 34881s 14 - 30854lV 3 - 1157621 3 s 12 - 2471991V 1 - 330131lV° - 280925lV 

- 1537051 7 s 8 - 830451 s s 7 - 1073851 9 s 6 - 14215U 10 s 5 - 1194591 n s 4 - 612421 12 s 3 

- 17534l 13 s 2 - 20481 14 s + 32l 15 )m 4 + 2sl(s + l) 2 (392s 12 + 29961s 11 + 9224lV° 

+ 145021 3 s 9 + 10803t 4 s 8 - 1352lV - 101461V - 70321 7 s 5 + 2203lV + 6782l 9 s 3 
+ 5064l 10 s 2 + 17561 n s + 232l 12 )m 2 — s 2 l 2 (s + l) 3 (s 2 + Is + l 2 ) 2 (20s 6 + 461s 5 

- 91lV - 238lV - 91lV + 461 5 s + 201 6 )] (AlOi) 
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/i 4 = -12 ml [201326592s(s - t)t{s + tf mf - 4194304st(s + t) 2 (291s 3 + 899ts 2 

+ 997t 2 s + 285 t 3 )mf + 524288(32s 8 + 3928ts 7 + 22769tV + 56744t 3 s 5 + 76938tV 
+ 59740t 5 s 3 + 25149t 6 s 2 + 4444t 7 s - 32 t 8 )m 2 c ° - 131072(224s 9 + 13750ts 8 
+ 84023tV + 231650t 3 s 6 + 369902t 4 s 5 + 371442t'V + 235061t 6 s 3 + 86390tV 
+ 13894t s s - 256t 9 )m 48 + 32768(672s 10 + 30316ts 9 + 191386tV + 564744tV 
+ 1001907t 4 s 6 + 1182124t 5 s 5 + 960430tV + 525736t 7 s 3 + 174749t 8 s 2 + 25168t 9 s 

- 896t 10 )m 46 - 8192(1120s n + 45670ts 10 + 298834tV + 931048t 3 s 8 + 1770223tV 
+ 2303172t 5 s 6 + 2179344t 6 s 5 + 1526640tV + 760913t 8 s 3 + 236598t 9 s 2 + 30486t 10 s 

- 1792t n )m 44 + 2048(1120s 12 + 48192ts n + 333351tV° + 1101380tV + 2213799tV 
+ 3049824t 5 s 7 + 3122706t 6 s 6 + 2517072tV + 1614829t 8 s 4 + 772192tV + 231987t 10 s 2 
+ 27212t n s - 2240t 12 )m 42 - 512(672s 13 + 35098ts 12 + 264391t 2 s n + 944150tV° 

+ 2026693t 4 s 9 + 2927376i 5 s 8 + 3093530tV + 2615244t 7 s 6 + 1916873t s s 5 + 1212706tV 
+ 586319t 10 s 3 + 174642t n s 2 + 19394t 12 s - 1792t 13 )m 40 + 128(224s 14 + 16740ts 13 
+ 142436t 2 s 12 + 565020t 3 s n + 1329468tV° + 2056868t 5 s 9 + 2236868 t 6 s 8 + 1858992tV 
+ 1367893t s s 6 + 1012724t 9 s 5 + 691232t 10 s 4 + 349104t 1 V + 103807t 12 s 2 + 11424t 13 s 

- 896t 14 )m 8 - 32(32s 15 + 4676ts 14 + 47458t 2 s 13 + 216720tV 2 + 577256t 4 s n 

+ 990974t 5 s 10 + 1149214tV + 936240t 7 s 8 + 606773tV + 438198t 9 s 6 + 396540t 10 s 5 
+ 310428t 1 V + 161311t 12 s 3 + 47012t 13 s 2 + 5248t 14 s - 256 t 15 )m 6 c + 8t(576s 15 
+ 8074ts 14 + 45396tV 3 + 142830t 3 s 12 + 281476t 4 s n + 360889tV° + 296852t 6 s 9 
+ 148525tV + 60588t s s 7 + 78887t 9 s 6 + 117204t 10 s 5 + 101911t n s 4 + 51660t 12 s 3 
+ 14268t 13 s 2 + 1600t 14 s - 32 t 15 )m 4 c - 2 st 2 (s + t) 2 (396s 12 + 2764ts n + 8634tV° 

+ 13982tV + 10851t 4 s 8 - 1158t 5 s 7 - 10180t 6 s 6 - 7280t 7 s 5 + 2283 t 8 s 4 + 6972t 9 s 3 
+ 5136t 10 s 2 + 1712t n s + 232 t 12 )m 2 + s 2 t 3 (s + t) 3 (s 2 + ts + t 2 ) 2 

(6s 6 + 17ts 5 - 105t 2 s 4 - 224t 3 s 3 - 62t 4 s 2 + 60t 5 s + 20t 6 ) 1 (AlOj) 
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h\ = -st (-4 m 2 + s + t) [622592m 34 (3s 2 + 2 st + 3t 2 ) - 4096m 32 (203s 3 - 12s 2 t 
+ 216st 2 + 203t 3 ) - 1024m 30 (420s 4 + 2357 sh + 2358s¥ + 1187st 3 + 258t 4 ) 

+ 256m 3 (2073s 5 + 8653s 4 t + 12881s¥ + 10805s¥ + 5215st 4 + 1263t 5 ) 

- 64 m 6 (3265s 6 + 14055s s t + 26073sV + 28978s¥ + 20175s¥ + 8319st 5 + 1807t 6 ) 

+ 48m 4 (774s 7 + 3736 sH + 8400s¥ + 11507s¥ + 10373s¥ + 6186s¥ + 2170st 6 

+ 396t 7 ) - 108m 2 (23s 8 + 142s 7 t + 392s¥ + 653s¥ + 742s¥ + 569s¥ + 296s¥ 

+ 88st 7 + lit 8 ) + 297st(s + t) (s 2 + st + t 2 ) 3 (Allb) 

h 2 = 8m 2 J4980736m 36 s - 16384m 34 (214s 2 + 13st + 271 2 ) + 4096m 32 (125s 3 

- 1442s 2 t - 799 st 2 + 216t 3 ) + 1024m 30 (1299s 4 + 7096s 3 t + 7922s¥ + 1845st 3 

- 756t 4 ) - 256m 8 (4029s 5 + 17654s 4 t + 25019s¥ + 14306s 2 * 3 + 1152st 4 - 15127 s ) 

+ 64m 6 (5497s 6 + 24150s 5 t + 40980s 4 * 2 + 33774s 3 * 3 + 11507s 2 * 4 - 1511st 5 - 18637 6 ) 

- 16m 4 (3672s 7 + 17571s 6 * + 35040s 5 * 2 + 36622s¥ + 19839s 3 * 4 + 2237s¥ - 3125st 6 

- 1404t 7 ) + 4 m 2 (918s 8 + 5400s 7 * + 13063s 6 * 2 + 16698s 5 * 3 + 11639s 4 * 4 + 2810s 3 * 5 

- 2283s¥ - 2057s* 7 - 594t 8 ) + 27*(-10s 8 - 40s 7 * - 63s 6 * 2 - 49s 5 * 3 + 2s 4 * 4 + 39s¥ 

+ 39s 2 * 6 + 18st 7 + 4* 8 ) (Allc) 

h :i = 8m 2 (4m 2 - s - *) [l024m 30 (27s 4 + 260s 3 * + 162s 2 * 2 + 260s* 3 + 27* 4 ) 

- 256m 8 (l62s 5 + 617s 4 * + 939s 3 * 2 + 1395s 2 * 3 + 617s* 4 + 162* 5 ) + 64m 6 (405s 6 
+ 1238s 5 * + 1841s 4 * 2 + 3670s 3 * 3 + 2717s¥ + 1172s* 5 + 405* 6 ) - 32m 4 (270s 7 

+ 642s 6 * + 1190s 5 * 2 + 2225s 4 * 3 + 2387s¥ + 1757s 2 * 5 + 723st 6 + 270* 7 ) + 4m 2 (378s 8 
+ 1193s 7 * + 2369s¥ + 3793s 5 * 3 + 4964s 4 * 4 + 4117s 3 * 5 + 2855s 2 * 6 + 1193s* 7 + 378t 8 ) 

- 27 (4s 9 + 18s 8 * + 39s 7 * 2 + 66s 6 * 3 + 83s¥ + 83s 4 * 5 + 66s¥ + 39s¥ + 18s* 8 + 4* 9 ) 

(Alld) 



K = 8m 2 [221184m 44 (s - t)(s + t) 2 - 2048m 42 (189s 4 - 110s 3 t - 162s¥ - 410st 3 

- 243t 4 ) + 512m 40 (567s 5 + 566s 4 t - 495s¥ - 2227s¥ - 2374st 4 - 945t 5 ) 

- 128m 3 (945s 6 + 2744s 5 t + 2585s 4 t 2 - 4606s¥ - 9379s¥ - 6778st 5 - 2079t 6 ) 

+ 32m 6 (918s 7 + 4957s 6 t + 8354s¥ - 709s¥ - 15709s¥ - 19272s¥ - 10811st 6 

- 2808t 7 ) - 8m 4 (486s 8 + 3942s 7 t + 9139s¥ + 4902s¥ - 11758s¥ - 23880s¥ 

- 21205s 2 t 6 - 9976st 7 - 2322t 8 ) + 2m 2 (108s 9 + 918s 8 t + 2906s¥ + 2147s¥ 

- 5453s¥ - 14923sY - 17387s¥ - 12136s¥ - 5032st 8 - 1080t 9 ) + 27t(2s 9 + 4s 8 t 

+ 13s¥ + 42s¥ + 83s¥ + 107s 4 t 5 + 92 s¥ + 53s 2 t 7 + 20st 8 + 4t 9 )] (Alle) 


g + g^ cc(¥j 81 ) + g: 


b"( 3 Pf ] ) = --- (Al2a) 

m b c s 2 t 2 (s — 4 m;?) 4 (t — Am 2 c ) ( s + ¥ (—4m;? + s + t) 2 


hi = -st (-4m 2 + s + 1) [l048576m 20 (51s 4 + 187s 3 t + 244s¥ + 127st 3 + 51t 4 ) 

- 262144m 48 (306s 5 + 1357s 4 t + 2375s¥ + 1905s¥ + 1007st 4 + 306t 5 ) 

+ 65536m 46 (911s 6 + 4526s 5 t + 9668s¥ + 10344s¥ + 7048s¥ + 3356st 5 + 831t 6 ) 

- 32768m 44 (807s 7 + 4441s 6 t + 10914s¥ + 14369s¥ + 12229s¥ + 7369s¥ 

+ 3126st 6 + 657t 7 ) + 8192 m 42 (888s 8 + 5519s 7 t + 15179s¥ + 23051s¥ + 23276s¥ 

+ 17136s¥ + 9319sY + 3604st 7 + 648t 8 ) - 2048m 40 (635s 9 + 4541s 8 t + 13692s¥ 

+ 22822s¥ + 25571s¥ + 21816s¥ + 14607s¥ + 7442s¥ + 2701st 8 + 405t 9 ) 

+ 256m 8 (627s 40 + 5 1 95s 9 t + 16915s¥ + 29976s 7 t 3 + 34868s¥ + 31408s 5 t 5 
+ 24488s¥ + 16106s 3 t 7 + 8015s 2 t 8 + 2755st 9 + 327t 40 ) - 64m 6 (210s n + 2085s 10 t 
+ 7585s¥ + 14784s 8 t 3 + 18290s¥ + 16790s¥ + 13650s 5 t 6 + 10810sV + 7164s 3 t 8 
+ 3345s 2 t 9 + 1005st 10 + 90t 44 ) + 16m 4 (35s 42 + 506s 44 t + 2338s 10 t 2 + 5745s¥ 

+ 9126s 8 t 4 + 10612s¥ + 9988s 6 t 6 + 8172s 5 t 7 + 5746s 4 t 8 + 3105s¥ + 1138sY° 

+ 246st 41 + 15t 42 ) — 4m 2 st (s 2 + st + t 2 ) 2 (46s 7 + 270s 6 t + 611s 5 t 2 + 761s¥ 

+ 681s¥ + 451s¥ + 170st 6 + 261 7 ) + lls 2 t 2 (s + t) 2 (s 2 + st + t 2 ) 4 (Al2b) 
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/i '2 = 32 m 2 4194304s(s — t) 2 tm 22 — 262144st (lls 3 — 44ts 2 + llt 2 s + 6 1 3 ) /n 20 

- 65536 (76s 6 + 237 ts 5 + 436tV - 153t 3 s 3 - 126tV + 118t 5 s + 36t 6 ) m 48 
+ 16384(416s 7 + 1681ts 6 + 2772 t 2 s 5 + 495t 3 s 4 - 892t 4 s 3 + 346tV + 998t 6 s 

+ 288 t ! )mf - 4096(940s 8 + 4547ts 7 + 8571t 2 s 6 + 4687t 3 s 5 - 470tV + 878t 5 s 3 
+ 4017t 6 s 2 + 3802t 7 s + 1008t s )m 34 + 1024(1120s 9 + 6286ts 8 + 13654tV + 11090t 3 s 6 
+ 3254tV + 4751t 5 s 4 + 10960t 6 s 3 + 12707tV + 8314t 8 s + 2016t 9 )m 42 - 256(740s 10 
+ 4700ts 9 + 11283t 2 s 8 + 11038tV + 5528t 4 s 6 + 9919fV + 21427t 6 s 4 + 26283 t 7 s 3 
+ 21360tV + 11286t 9 s + 2520t 10 )m 30 + 64(256s n + 1780ts 10 + 4120tV + 3090t 3 s 8 
+ 544t 4 s 7 + 7817t 5 s 6 + 24064t 6 s 5 + 34463tV + 31912t 8 s 3 + 21376tV + 9738t 10 s 
+ 2016t 11 )m 8 - 16(36s 12 + 244ts n - 40tV° - 2560t 3 s 9 - 5744t 4 s 8 - 1877tV 
+ 11279t 6 s 6 + 24269t 7 s 5 + 27618tV + 21672t 9 s 3 + 12769t 10 s 2 + 5218t 11 s + 1008t 12 )m 6 
+ 8t(—6s 12 - 192ts n - 1024t 2 s 10 - 2580t 3 s 9 - 3274t 4 s 8 - 1283t 5 s 7 + 2616t e s 6 
+ 5395tV + 5517t 8 s 4 + 3840tV + 2057t 10 s 2 + 786 t u s + 144t 12 )m 4 + t 2 (52s 12 
+ 396ts n + 1396t 2 s 10 + 2700t 3 s 9 + 2943t 4 s 8 + 1412t 5 s 7 - 850t 6 s 6 - 1920tV - 1653t 8 s 4 

- 916tV - 452t 10 s 2 - 188t n s - 36 t 12 )m 2 c - st 3 (s 2 + ts + t 2 ) 2 (s 7 + 4 ts 6 - 4 t 2 s 5 

- 26t 3 s 4 - 22t 4 s 3 + 4tV + 9t 6 s + 2t 7 )} (Al2c) 
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h 3 = 32 m 2 c [262144m I s st (28s 4 + 93 s 3 t + 126s¥ + 73 st 3 + 28t 4 ) - 16384m 46 st(653s 5 
+ 2802s 4 t + 4924s¥ + 4254s¥ + 2407st 4 + 688t 5 ) - 4096m 44 (116s 8 - 1263s 7 t 

- 8068s 6 t 2 - 18556s 5 * 3 - 21164s 4 * 4 - 15711s 3 * 5 - 7768s¥ - 1718st 7 + 36t 8 ) 

+ 1024m 42 (656s 9 + 369s 8 t - 9566s 7 * 2 - 32833s 6 * 3 - 49800s 5 * 4 - 46665s 4 * 5 

- 30158s 3 * 6 - 12591s 2 * 7 - 1916s* 8 + 216* 9 ) - 256m 40 (1540s 10 + 5543s 9 * + 1997s 8 * 2 

- 22185s 7 * 3 - 53859s 6 * 4 - 66855s¥ - 54509s 4 * 6 - 29995s 3 * 7 - 9693s 2 * 8 - 92s* 9 

+ 540* 10 ) + 64m 8 (1920s n + 9816s 10 * + 18914s 9 * 2 + 13315s 8 * 3 - 10042s 7 * 4 - 33052s 6 * 5 

- 39132s 5 * 6 - 27447s 4 * 7 - 11120s 3 * 8 - 576s 2 * 9 + 2236s* 10 + 720* 11 ) - 16m 6 (1340s 12 
+ 8240s 11 * + 21729s 10 * 2 + 32953s 9 * 3 + 32174s 8 * 4 + 21270s 7 * 5 + 9410s 6 * 6 + 3455s 5 * 7 
+ 3594s 4 * 8 + 4948s 3 * 9 + 4629s 2 * 10 + 2550s* 11 + 540* 12 ) + 8m 4 (248s 13 + 1724s 12 * 

+ 5352s 11 * 2 + 10237s 10 * 3 + 14018s 9 * 4 + 15183s 8 * 5 + 13805s 7 * 6 + 10710s 6 * 7 + 7198s 5 * 8 
+ 4463s 4 * 9 + 2707s 3 * 10 + 1507s 2 * 11 + 604s* 12 + 108* 13 ) - m 2 (s + *) 2 (76s 12 + 412s 11 * 

+ 1004s 10 * 2 + 1740s 9 * 3 + 2541s 8 * 4 + 3200s 7 * 5 + 3202s 6 * 6 + 2320s¥ + 1141s¥ 

+ 440s¥ + 204s 2 * 10 + 132s* 11 + 36* 12 ) — s*(s + *) 3 (s 2 + st + * 2 ) 2 (2s 6 + 5s' 5 * — 8s¥ 

- 23s¥ - 8s¥ + 5s* 5 + 2* 6 )] (Al2d) 
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h A = 32m 2 [65536st (59s 4 + 150ts 3 + 252 t 2 s 2 + 182t 3 s + 53t 4 ) m 18 - 16384(96s 7 
+ 690ts 6 + 1748t 2 s 5 + 2626 t 3 s 4 + 2313tV + 1154fV + 213t 6 s - 36t 7 )m 46 
+ 4096(536s 8 + 3159ts 7 + 8087t 2 s® + 12811t 3 s 5 + 12362tV + 7490t 5 s 3 + 2561t 6 s 2 

- 26 t 7 s - 252t 8 )m 44 - 1024(1240s 9 + 7337 ts 8 + 19822tV + 33719t s s 6 + 36841t 4 s 5 
+ 25464t 5 s 4 + 10399t 6 s 3 + 1164tV - 1706t 8 s - 756t 9 )m 42 + 256(1520s 10 + 9504ts 9 
+ 27735t 2 s 8 + 50970t 3 s 7 + 62054t 4 s 6 + 49255fV + 22995t®s 4 + 2530t 7 s 3 - 4972t 8 s 2 

- 4387t 9 s - 1260t 10 )m 40 - 64(1040s n + 6961ts 10 + 22156tV + 44286t 3 s 8 + 58938t 4 s 7 
+ 51617t 5 s 6 + 26090t 6 s 5 + 904t 7 s 4 - 10831t s s 3 - 10450t 9 s 2 - 5487t 10 s - 1260t n )m 8 
+ 16(376s 12 + 2682ts n + 9390t 2 s 10 + 20757t 3 s 9 + 30044t 4 s 8 + 27253t'V + 11565t 6 s 6 

- 5978tV - 15058t 8 s 4 - 14380t 9 s 3 - 9099t 10 s 2 - 3812t n s - 756t 12 )m 9 - 4(56s 13 
+ 408ts 12 + 1632t 2 s n + 4317t 3 s 10 + 7054t 4 s 9 + 5921t 5 s 8 - 796tV - 8944t 7 s 6 

- 12637t 8 s 5 - 11090t 9 s 4 - 7241t 10 s 3 - 3756t n s 2 - 1396t 12 s - 252 t 13 )m 4 c - t(8s 13 

- 186t 2 s n - 472t 3 s 10 - 9 t 4 s 9 + 1935 t 5 s 8 + 4462t®s 7 + 5508t 7 s 6 + 4381t s s 5 + 2457tV 
+ 1116t 10 s 3 + 504t n s 2 + 196t 12 s + 36 t l3 )m 2 — st 2 (s 3 + 2 ts 2 + 2 t 2 s + t 3 )~ 

(s 6 + 3 ts 5 - 9 t 2 s 4 - 22 t 3 s 3 - 6 t 4 s 2 + 6 t 5 s + 21 6 ) 1 (Al2e) 

9 + 9 -> cc( 3 Sf ] , 3 D [8] ) + g: 

b" ( 3 4 81 , 3 ^i 81 ) =- -i= -4-- (Al3a) 

54-\/T5 mist (s — 4 m;?) 4 (t — 4 m%) (s + t) 4 (—4 m 2 c + s + t) 

hi = st (s — 4m 2 ) (4m 2 — t) (s + t ) (4m 2 — s — t) 49152m 46 (1261s + 1459t) 

- 1024m 44 (34157s 2 + 40400st + 42419t 2 ) + 512m 42 (5978s 3 - 29869s 2 t - 21472st 2 
+ 11999t 3 ) + 64m 40 (10841s 4 + 214882s 3 t + 450086s 2 t 2 + 196198st 3 + 5225t 4 ) 

+ 32m 8 (19811s 5 - 64363s 4 t - 264691s 3 t 2 - 265735s 2 * 3 - 67162s* 4 + 19244* 5 ) 

- 4m® (86162s 6 + 168328s 5 * + 43407s 4 * 2 - 172806s 3 * 3 + 34263sV + 163270s* 5 
+ 85568* 6 ) + m 4 (57024s 7 + 220986s 6 * + 406490s 5 * 2 + 459356s 4 * 3 + 459698s 3 * 4 
+ 407102s 2 * 5 + 221256s* 6 + 57024* 7 ) - 162m 2 (22s 8 + 128s 7 * + 331s 6 * 2 + 535s 5 * 3 

+ 632s 4 * 4 + 535s 3 * 5 + 331s 2 * 6 + 128s* 7 + 22t 8 ) + 405s*(s + *) (s 2 + st + * 2 ) 3 1 (Al3b) 
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h 2 = -2 m\t (4m 2 - t) [l31072m 48 s(19027s + 38911*) - 32768m 46 (34301s s 

+ 175284s 2 * + 117007s* 2 + 1620* 3 ) - 4096m 44 (81504s 4 - 265155s 3 * - 587452s 2 * 2 

- 154537s* 3 - 25920* 4 ) + 1024m 42 (340730s 5 + 674275s 4 * + 265516s 3 * 2 + 143007s 2 * 3 
+ 77548s* 4 - 90720* 5 ) - 256m 40 (443414s 6 + 1706395s 5 * + 2270361s 4 * 2 + 1548543s 3 * 3 
+ 710607s 2 * 4 - 161108s* 5 - 181440* 6 ) + 64m 8 (367548s 7 + 1961043s 6 * + 3676724s 5 * 2 
+ 3126628s 4 * 3 + 1228754s 3 * 4 - 268787s 2 * 5 - 686126s* 6 - 223560* 7 ) - 16m 6 (235756s 8 
+ 1490664s 7 * + 3518089s 6 * 2 + 3849434s 5 * 3 + 1701248s 4 * 4 - 421485s 3 * 5 - 1160629s 2 * 6 

- 731477s* 7 - 168480* 8 ) + 4m 4 (103680s 9 + 755956s 8 * + 2167956s 7 * 2 + 3098173s 6 * 3 
+ 2077386s 5 * 4 + 205573s 4 * 5 - 775916s 3 * 6 - 761172s 2 * 7 - 345856s* 8 - 71280* 9 ) 

- m 2 (20736s 40 + 207360s 9 * + 737224s 8 * 2 + 1361080s 7 * 3 + 1370281s 6 * 4 + 690475s 5 * 5 
+ 23643s 4 * 6 - 195151s 3 * 7 - 151624s 2 * 8 - 60912s* 9 - 12960* 40 ) + 162s*(32s 9 + 160s 8 * 
+ 359s 7 * 2 + 476s 6 * 3 + 409s 5 * 4 + 233s 4 * 5 + 96s 3 * 6 + 39s 2 * 7 + 16s* 8 + 4* 9 )] (Al3c) 

h 3 = —2m 2 (s + *) (4m 2 - s - *) [65536m 46 s*(6701s + 6899*) + 8192m 44 (1620s 4 

- 63459s 3 * - 138920s 2 * 2 - 70317s* 3 + 1620* 4 ) - 1024m 42 (22680s 5 - 201895s 4 * 

- 495519s 3 * 2 - 539529s 2 * 3 - 244321s* 4 + 22680* 5 ) + 256m 40 (68040s 6 + 28875s 5 * 

- 140260s 4 * 2 - 209950s 3 * 3 - 209740s 2 * 4 - 717s* 5 + 68040* 6 ) - 64m 8 (113400s 7 
+ 353683s 6 * + 595357s 5 * 2 + 986984s 4 * 3 + 1022588s 3 * 4 + 637657s 2 * 5 + 359587s* 6 

+ 113400* 7 ) + 16m 6 (110160s 8 + 383177s 7 * + 750828s 6 * 2 + 1345133s 5 * 3 + 1831816s 4 * 4 
+ 1470557s 3 * 5 + 821388s 2 * 6 + 389549s* 7 + 110160* 8 ) - 4m 4 (58320s 9 + 205510s 8 * 

+ 354804s 7 * 2 + 556953s 6 * 3 + 851841s 5 * 4 + 878157s 4 * 5 + 598875s 3 * 6 + 370788s 2 * 7 
+ 205888s* 8 + 58320* 9 ) + m 2 (s + *) 2 (12960s 8 + 16848s 7 * + 262s 6 * 2 - 23497s 5 * 3 
+ 8468s 4 * 4 - 26953s 3 * 5 - 8s 2 * 6 + 16848s* 7 + 12960* 8 ) + 162s*(4s 9 + 24s 8 * + 75s 7 * 2 
+ 120s 6 * 3 + 149s 5 * 4 + 149s 4 * 5 + 120s 3 * 6 + 75s 2 * 7 + 24s* 8 + 4* 9 ) (Al3d) 
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h 4 = m 2 c [3145728m 20 st(2357s - 3083t) - 131072m 48 (1620s 4 + 105283s 3 t + 8642 s 2 t 2 

- 96641st 3 - 1620t 4 ) + 16384m 46 (25920s 5 + 509621s 4 t + 896167sV - 85037s¥ 

- 245935st 4 - 32400t 5 ) - 8192m 44 (45360s 6 + 320292 s b t + 1056932sV + 856843s¥ 

+ 146706s¥ + 583 st 5 - 712801 6 ) + 1024m 42 (181440s 7 + 599428s 6 t + 1988277s'V 

+ 2840253s¥ + 1083297sV + 141039s¥ - 376166st e - 362880t 7 ) 

- 256m 40 (223560s 8 + 623352s 7 t + 977271s¥ + 992272s¥ - 1796934s¥ 

- 3815074s¥ - 2849597s¥ - 1832690st 7 - 586440t 8 ) + 64m 8 (168480s 9 + 592859s s t 
+ 679396s¥ - 438617s¥ - 4960050s¥ - 10323584s¥ - 9863307s¥ - 6026378sV 

- 2744639st 8 - 615600t 9 ) - 16m 6 (71280s 40 + 312714s 9 t + 668005s¥ + 413842s¥ 

- 2602893s¥ - 8848230s¥ - 12123493s¥ - 9435072s¥ - 5216495s¥ - 2055090st 9 

- 408240t 10 ) + 4m 4 (12960s n + 49248s 10 t + 294940sV + 872129s¥ + 834712s¥ 

- 1612967s¥ - 5014855s¥ - 5770958sV - 4063097s 3 * 8 - 2154136s 2 * 9 - 805760st 40 

- 155520* 44 ) + m 2 t(7776s 44 - 55728s 40 * - 413458s 9 * 2 - 1003935s 8 * 3 - 1092466s 7 * 4 

- 310362s 6 * 5 + 540758s 5 * 6 + 744113s 4 * 7 + 542710s 3 * 8 + 309728s 2 * 9 + 124416s* 10 

+ 25920* 14 ) + 324s* 2 (s + t) 2 (8s 8 + 44s 7 * + 83 s 6 t 2 + 85 s 5 t 3 + 64s 4 * 4 + 37 s 3 t 5 + 31s 2 * 6 
+ 16s* 7 + 4t 8 )j (Al3e) 

q(q) +9 cc( 3 Pj 1] ) +q(q): 


a <4<?W3p[ih _ 32¥a 3 

1 4 ; 81m3(s + t) 4 

(Al4a) 

gi — — (s + *) [-4m 2 (s + t) + s 2 + t 2 ] 

(Al4b) 

g -2 = —16 m 2 s 

(Al4c) 

33 = 0 

(Al4d) 

g 4 = -8m 2 (s -*) 

(Al4e) 

uq(q)a(3pW\ _ ■*'® 7r “s 

4 405m® (4m;? — s) (s + t) 5 

(Al4f) 
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h\ = — (s +1) 16m 4 (s 2 + 22st + 21t 2 ) — 4m 2 (12s 3 + 43s 2 t — 8st 2 + 21t 3 ) 

+ 11a (s 3 + s 2 t + st 2 + t 3 ) (Al4g) 

h 2 = 16m 2 s [4m 2 (s + 41 1) — s(s + t)] (Al4h) 


ti, = 0 


(Al4i) 


hi = 8m 2 [4m 2 (s 2 + 50st — 31t 2 ) — s 3 + st 2 


?(?) + 3 ->• cc( 3 Pj 1] ) + g(g): 


a ?(9)fl(3p[l]^ 


32vr 3 a 3 

1215mj?(s + t ) 4 (—4m;? + s + t) 


(A14j) 


(Al5a) 


g\ = 16 m 4 (l9s 2 + 30 st + 19t 2 ) — 8m 2 (s 3 + 16s 2 t + 16st 2 + t 3 ) + (s + t) 2 (s 2 + t 2 ) 


(Al5b) 


r/ 2 = 48m 2 (32 m 4 + 4 m 2 (s + 2t) — s 2 — st + 2t 2 ) 


(Al5c) 


93 = 96m 2 (s + t) 2 


(Al5d) 


,r/ 4 = —24m 2 (4m 2 (s — t) + s 2 — 4st — 51 2 ) 


(Al5e) 


l,q(<i)9f3p[i]\ _ __ 

6075m® (4m;? — s) (s + t) 5 (—4m;? + s + t) 2 


(Al5f) 


ti, = - (4m 2 - s - t) [320m® (l9s 3 + 141s 2 t + 153st 2 + 95t 3 ) 

+ 16 m 4 (73s 4 + 39s 3 t - 923s 2 t 2 - 931st 3 - 42t 4 ) 

+ 4m 2 (s + t) 2 (l9s 3 + 139s 2 t — st 2 + 17t 3 ) — 7s(s + t) 3 (s 2 + t 2 ) 


(Al5g) 


= 48m 2 2560m 3 (s — 71) — 64m® (49s 2 — 3st + 4t 2 ) 

+ 16m 4 (37s 3 + 113s 2 t + 180st 2 + 8t 3 ) + 4m 2 (5s 4 - 21s 3 t - 69s 2 t 2 + 23st 3 + 66t 4 ) 

— s(s + t) 2 (3s 2 — st + 6t 2 ) (Al5h) 
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h 3 = —96m 2 (s + t) 2 80m 4 (3s + 7t) — 4m 2 (16s 2 + 47st + 31t 2 ) + s(s + t) 2 (Al5i) 


h A = 24m : 


320 ml (s 2 + 10 st - 7 t 2 ) - 80m 4 (s 3 + 8 sH + 31st 2 + 24t 3 ) 
— 4m 2 (s + t) 2 (s 2 + 26st — 155t 2 ) + s(s — 5t)(s + t) 3 

q(q) + g ->• cc( 3 5f ] ) + g(g): 

6?(5)fl(3g-| 8 ]) = 


7T 3 a 3 


162m^st (s — Ami) ( s + ^) 3 (—4mj? + s + t) 


(Al5j) 


(Al6a) 


h = 


- (4m 2 - s - t) 128ml (20s 3 + 69s 2 t - 39st 2 + 20t 3 ) - 32m 4 (40s 4 + 104s 3 t 
+ 45s 2 t 2 + st 3 + 20t 4 ) + 4m 2 (l08s 5 + 175s 4 t + 77 s 3 t 2 + 207s 2 t 3 + st 4 + 20t 5 ) 


11s (4s 5 + 3s 4 t + 7 s 3 t 2 + 7s 2 t 3 + 3st 4 + At 5 ) 


(Al6b) 


ho = —16m 2 


576m 5 (s - t) 2 + 16m 4 (26s 3 + 105s 2 t - 3 st 2 + 26t 3 ) - 4m 2 (55s 4 


+ 155s 3 t + 99s 2 t 2 + 34st 3 + 35t 4 ) + s (20s 4 + 44s 3 t + 21s 2 t 2 + 44st 3 + llt 4 )(AT6c) 


h 3 = — 32m 2 (4m 2 — s — t) 4m 2 (44s 3 + 87s 2 t — 21st 2 + 44t 3 ) 
— 5s (4s 3 + 3s 2 t + 3st 2 + 4t 3 ) 


(Al6d) 


hd = —8m 2 


16m 4 (97s 3 + 147s 2 t - 15st 2 + 791 3 ) - 4m 2 (146s 4 + 265s 3 t + 189s 2 t 2 


+ 77st 3 + 79t 4 ) + s (49s 4 + 70s 3 t + 60s 2 t 2 + 70st 3 + 31t 4 ) 
q(q)+9~> cc( 3 Pf ] ) + q(q): 

M(q)9(3pM\ _ ^ JQs M(gbf3p[8]\ 
v ; 128a V J ’ 

The coefficients hf, are the same as those in Eq. (1A6jl . 
q(q) + 9 -t cc( 3 Sj 8] , 3 A>! 81 ) + q{q): 


(Al6e) 


(A17) 


jjIsisjasIsl SpIBl) _ 


7r 3 a 3 


27v / 15?« 5 st(s + t) 4 (—Ami + s + t) 


(Al8a) 


hi = — (s + t) (—4m 2 + s + t) 160m 4 (4s 3 + 3s 2 t + 3st 2 + 4t 3 ) — 4m 2 (40s 4 + 97s 3 t 


+ 6s 2 t 2 + 97st 3 + 40t 4 ) + 5 (4s 5 + 3s 4 t + 7sY + 7s 2 t 3 + 3st 4 + 4t 5 ) 


(Al8b) 
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h 2 = 8 m 2 3456m® (s — t ) 2 — 144 m 4 (4s 3 — s 2 t — 16st 2 + 13t 3 ) + 4m 2 (22s 4 + 25s 3 t 
+ 114s 2 t 2 + 7 st 3 + 112t 4 ) - 40s 5 - 128s 4 t - 157sV - 103s 2 t 3 - 83st 4 - 49tfl\18c) 


h :i = 8m 2 (s + t ) 2 (l07s 2 — 74st + 107t 2 ) (4m 2 — s — t) (Al8d) 

h 4 = 8 m 2 c [216m 4 (3s 3 - 5s 2 t + 5 st 2 - 3t 3 ) + 2m 2 (17s 4 + 158s 3 t + 66s 2 t 2 


+ 122st 3 + 197t 4 ) - (s + t ) 2 (49s 3 + 48s 2 t - 15st 2 + 58t 3 ) 

(Al8e) 

q + q —* cc( 3 P 1 [1] ) + g: 


a ™( 3 P [1] ) - 2 56vr 3 a 3 

1 243m (s — 4ml) 4 

(Al9a) 

gi — — (4m 2 — s) (4m 2 (s + 2t) — s 2 — 2 st — 2 1 2 ) 

(Al9b) 

g 2 = 16m 2 (4m 2 — s — t) 

(Al9c) 

g 3 = 16 m 2 c t 

(Al9d) 

r/ 4 = 32 m 4 — 8m 2 s 

(Al9e) 

jjW^3p[i]\ _ 128vr 

1 1215m® (4m;? — s) 5 

(Al9f) 

h 1 = (4m 2 - s) [l6m 4 (31s + 42t) - 8m 2 (21s 2 + 32 st + 21 1 2 ) + 11s (s 2 

+ 2 st + 2 1 2 ) ] 


(Al9g) 

h 2 = —16m 2 (124m 2 — s) (4m 2 — s — t) 

(Al9h) 

h 3 = 16 m 2 c t (s — 124m 2 ) 

(Al9i) 

h 4 = —8m 2 (4m 2 — s) (l24m 2 — s) 

(Al9j) 
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q + q —> cc( 3 P^) + g: 


= — 2^1 — 
3645mj!s (s — 4 m%) 


g x = -4608m® + 2304m®(s + t)~ 16m 3 (19s 2 + 66st + 36t 2 ) 
+ 8m 2 s (s 2 + 16si + 15t 2 ) — s 2 (s 2 + 2st + 2t 2 ) 


(A20a) 


(A20b) 


g -2 = —48m 2 (32m 3 + 4m 2 (s + 2t) — s 2 — st + 2t 2 ) 


(A20c) 


g 3 = —48m 2 (96 m 3 — 24m 2 (s + t) + 2s 2 + 5 st + 2t 2 ) 


(A20d) 


= —24m 2 (96 m 3 — 4m 2 (s + At) — s 2 + Ast + 4t 2 ) 


6««( 3 P 2 [1] ) =- 128?r3Q ^ - F 

18225m 3 s 2 (4m 2 - s) 5 


= s [276480m 30 - 512m®(221s + 246t) + 64m® (93s 2 + 950st + 492t 2 ) 
+ 16m 3 s (53s 2 - 598 st - 458t 2 ) + 4m 2 s 2 (33s 2 + 154st + 140t 2 ) 

- 7s 3 (s 2 + 2st + 2t 2 ) 

h 2 = 48m 2 [2048m® + 128m®(19s - 8 1) + 16m 3 (57s 2 + 54st + 8 t 2 ) 

— 4m 2 s (20s 2 + st — 58t 2 ) — s 2 (3s 2 + 7 st + 10t 2 ) 

h 3 = 48m 2 [9600m®s - 32m 3 (58s 2 + 77 st - 4t 2 ) 

+ 4m 2 s (52s 2 + 137st + 58t 2 ) - s 2 (6s 2 + 13st + 10t 2 ) 

h A = 24m 2 128m®(85s - St) + 16m 3 (3s 2 - lOOst + 16t 2 ) 

- 8m 2 s (3s 2 - 68 st - 58 t 2 ) - s 2 (7s 2 + 20st + 20t 2 ) 


(A20e) 


(A20f) 


(A20g) 


(A20h) 


(A20i) 


(A20j) 


q + q ->• cc( 3 A( 81 ) + g: 


b ?q ( 3 sP) = 


47r 3 a 3 


2A3stml (s — Ami) {—Ami + s + t) 


(A21a) 
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hi — s [64m 4 — 4 m 2 (8s + 9 1) + 4s 2 + 9 st + 91 2 ] 704m® + 16m 4 (s — 221) 

+ 4m[ (23s 2 + 44st + 22t 2 ) - 11s (s 2 + 2st + 2 1 2 ) 1 (A21b) 


/i 2 = 16m 2 [64m®(77s - 18t) - 96m* (29s 2 + 24st - 12t 2 ) 

+ 36m 2 (l3s 3 + 22 s 2 t + 9st 2 - 10t 3 ) - 20s 4 - 36s 3 t - 9sV + 54st 3 + 36t 4 (A21c) 

h 3 = 16m 2 4352m®s — 48m 4 (49s 2 + 48st — 6t 2 ) 

+ 36m 2 (10s 3 + 16s 2 t + 9st 2 - 6t 3 ) - 11s 4 + 45s 2 t 2 + 90st 3 + 36t 4 ] (A21d) 



(A21e) 


(A22a) 


hi = s (4m 2 — s) 3264m® + 16m 4 (7s — 30 1) 

+ 4m 2 (33s 2 + 52 st + 30t 2 ) - 11s (s 2 + 2 st + 21 2 ) 1 (A22b) 


h 2 = 32m 2 [ll52m® + 48m 4 (5s - 12 1) + 8m 2 t(5s + 9 t) + s (s 2 + 2 st + 4 1 2 ) 1 (A22c) 

/i 3 = 32m 2 464m 4 s + 8m 2 (s 2 + 9 st + 91 2 ) + s (3s 2 + 6 st + 4t 2 ) (A22d) 

h A = 32m 2 16m 4 (19s — 18 1) + 8m 2 (7s 2 + 7 st + 9t 2 ) + s(s + 2t) 2 
q + g —> cc( 3 Af', 3 Zl[ 8 ') + 5 c 

6 *^, =-=-- 

81\/T5 m b c st (s — 4ml) (—4ml + s + t) 


(A22e) 


(A23a) 


hi = s (s - 4m 2 ) [20480m) 0 - 1024m®(15s + 28 t) + 64m®(80s 2 + 331st 

+ 310t 2 ) - 16m 4 (80s 3 + 411s 2 t + 705st 2 + 396t 3 ) + 4m 2 (60s 4 + 277s 3 t + 570s 2 t 2 
+ 576st 3 + 198t 4 ) - 5s (4s 4 + 17s 3 t + 35s 2 t 2 + 36st 3 + 18t 4 ) 1 (A23b) 
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h 2 = -8m 2 512m 3 (56s - 45 1) - 64 m\ (322s 2 + 117 st - 360t 2 ) + 48m 4 (116s 3 
+ 120s 2 t - 75 st 2 - 150t 3 ) + m 2 (-712s 4 - 900s 3 t + 144s 2 t 2 + 936st 3 + 720t 4 ) 

+ 40s 5 + 72 s 4 t + 45 s 3 t 2 + 36st 4 ] (A23c) 


h-t = —8m 2 


10240m®s - 64m 5 (115s 2 + 99st - 90t 2 ) + 48m 4 (53s 3 + 48s 2 t - 39st 2 


- 90t 3 ) - 4m 2 (133s 4 + 207s 3 t + 18s 2 t 2 - 342st 3 - 180t 4 ) 
+ s (49s 4 + 162s 3 t + 261s 2 t 2 + 144st 3 + 367 4 ) 


(A23d) 


h 4 = —72m 2 (—4m 2 + s + 2t)“ 40?n 4 (s — 2 1) + 2m 2 (9s 2 + 8 st + 10t 2 ) 
s (—s 2 4~ st 4" a) 


(A23e) 
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